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ABSTRACT 


In the formulations of earlier Displacement Transfer Functions for structure shape predictions, the 
surface strain distributions, along a strain-sensing line, were represented with piecewise linear functions. 
To improve the shape-prediction accuracies, Improved Displacement Transfer Functions were formulated 
using piecewise nonlinear strain representations. Through discretization of an embedded beam 
(depth-wise cross section of a structure along a strain-sensing line) into multiple small domains, 
piecewise nonlinear functions were used to describe the surface strain distributions along the discretized 
embedded beam. Such piecewise approach enabled the piecewise integrations of the embedded beam 
curvature equations to yield slope and deflection equations in recursive forms. The resulting Improved 
Displacement Transfer Functions, written in summation forms, were expressed in terms of beam 
geometrical parameters and surface strains along the strain-sensing line. By feeding the surface strains 
into the Improved Displacement Transfer Functions, structural deflections could be calculated at multiple 
points for mapping out the overall structural deformed shapes for visual display. The shape-prediction 
accuracies of the Improved Displacement Transfer Functions were then examined in view of 
finite-element-calculated deflections using different tapered cantilever tubular beams. It was found that by 
using the piecewise nonlinear strain representations, the shape-prediction accuracies could be greatly 
improved, especially for highly-tapered cantilever tubular beams. 

NOMENCLATURE 

A = x - 4 s l + s i+l )/2A l , 1/in. 

B = (£f_i - + £ i+1 )/2(A/) 2 , 1/in 2 

C depth factor (vertical distance from neutral surface to outermost fiber of bottom surface of 
uniform beam), in. 

c i = c(x z ), depth factor of nonuniform beam at strain-sensing station, X-, in. 
c n value of c i at free end (beam tip), x = x n = / , in. 
c Jc 0 depth ratio, no dimension 

c(x) depth factor (vertical distance from neutral surface to bottom surface outermost fiber of 
nonuniform beam) at axial location, X, in. 
c 0 value of c i at fixed end (beam root), x = x 0 = 0 , in. 

C = (c- - c i X ) A/ , in/in 

CG center of gravity 

di chord- wise distance between front and rear strain-sensing stations at { x i? x' }, in. 

D = [(c, / c i X ) — 1]/A /, 1/in. 

DLL design limit load 

i =0,1,2,3 , strain-sensing station identification number 

j dummy index 

/ length of beam, in. 

n index for the last span-wise strain-sensing station (or strain-sensing station density) 

P force, lb 

RRF rotated reference frame, aligns with the wing’s out-of-plane direction 
R(x) radius of curvature of beam elastic curve, in. 

SPAR Structural Performance And Resizing (finite-element computer program) 
t tubular beam wall thickness, in. 

x, y Cartesian coordinates (x in beam axial direction, y in lateral direction), in. 
x t axial coordinate or symbol of /- th strain-sensing station, in. 
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y(x) beam deflection in y-direction at axial location, X , in. 

y. = _y(x-) , beam deflection in y-direction at axial location, x = x n in. 

y^ PAR beam-tip deflection calculated from SPAR, in. 

A / = (x- — X £ l ) = / / n, distance between two adjacent strain-sensing stations { x i l? x-}, in. 

A y„ =\y n -yT l beam-tip shape-prediction error (difference between predicted and 

SPAR-generated beam-tip deflections), in. 
a =tan [(c 0 -c n )/l], beam taper angle, deg 

s(x) surface bending strain at axial location, X, in/in 
S t =s{x,), surface bending strain at strain-sensing station, x-, in/in 
0{x) slope of deformed beam at axial location, X, rad or deg 

o, =0(x,), slope of deformed beam at axial location, x = x /? rad or deg 
^ = (x — X- 1 ) , local axial coordinate with origin at x M , in. 

<j> i cross-sectional twist angle, rad or deg 
( )' quantities associated with rear strain-sensing line 


INTRODUCTION 


In the formulations of earlier Displacement Transfer Functions (refs. 1-5) for structure shape 
predictions using distributed surface strains, the embedded beam (depth-wise cross-section of the 
structure along a strain-sensing line) was first discretized into multiple small domains of equal length. 
Through discretization, the axial distributions of both beam depth and surface strain within each small 
domain could be represented with piecewise linear functions. The discretization approach enabled the 
piecewise integrations of the beam curvature equation (second-order differential equation) over the small 
domains in closed forms to yield beam slope and deflection equations in recursive forms. Each set of 
recursive slope and deflection equations was then combined into a single dual summation-form deflection 
equation called a Displacement Transfer Function (refs. 1-5). The Displacement Transfer Function was 
then written in terms of geometrical parameters of discretized embedded beam and surface strains 
obtained at strain-sensing stations (at domain junctures) evenly spaced along the strain-sensing line. By 
inputting the surface strains from multiple strain-sensing lines, the associated Displacement Transfer 
Functions can convert the surface strains into deflections at multiple points, enabling the mapping of 
overall structural deformed shapes for visual display. Thus, the Displacement Transfer Functions with an 
accompanying surface strain-sensing system [for example, fiber optic strain sensors (refs. 6-8)], has 
created a new structure shape-sensing technology, called Method for Real-Time Structure Shape-Sensing 
(U.S. Patent Number 7,520,176) (ref. 9), which could be very useful for in-flight deformed shape 
monitoring of flexible wings and tails, such as those often employed on unmanned flight vehicles by the 
ground-based pilot for maintaining safe flights. In addition, the wing shape being monitored in real time 
could then be input to the aircraft control system for aeroelastic wing shape control. 

Similar to the structure deformed shape predictions, the distributed surface strains can also be fed into 
the special stiffness and Load Transfer Functions described in Process for Using Surface Strain 
Measurements to Obtain Operational Loads for Complex Structures (U.S. Patent No. 7,715,994) (ref. 10), 
for predicting the operational loads on a flight vehicle. 
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The previously-developed Displacement Transfer Functions were successfully validated for their 
accuracies in view of finite-element analysis of different sample structures such as cantilever tubular 
beams (uniform, tapered, slightly tapered, and step-wisely tapered); two-point supported tapered tubular 
beams, flat panels, and tapered wing boxes (un-swept and swept) (refs. 1, 2, 4, and 5); and the 
doubly-tapered wing of the NASA Ikhana Unmanned Science and Research Aircraft System (a modified 
Predator-B unmanned aerial system, General Atomics Aeronautical Systems Inc., San Diego, California), 
(ref. 3). 

The earlier Displacement Transfer Functions were formulated for straight embedded beams. With the 
introduction of empirical curvature-effect correction terms, however, those Displacement Transfer 
Functions could be used for the shape predictions of slender curved structures with different arc-angles up 
to 360 deg (a complete circle) (ref. 11). 

All the earlier Displacement Transfer Functions were formulated using piecewise linear 
representations of the actual strain distributions. For the shape predictions of tapered beam structures, for 
which the strain gradients change rapidly (highly nonlinear strain distributions) (ref. 3), the use of 
piecewise linear strain representations could result in a loss of accuracy. One way to avoid this loss of 
shape-prediction accuracy is to increase the domain density (that is, use increasing strain-sensing stations) 
so that the piecewise linear strain representation can maintain a good approximation. An alternative 
approach is to introduce piecewise nonlinear representations of the actual strain distributions. 

This technical paper discusses the mathematical formulations of the Improved Displacement Transfer 
Functions using the piecewise nonlinear strain representations. A family of tapered cantilever tubular 
beams with different taper angles was used to study the shape-prediction accuracies of the Improved 
Displacement Transfer Functions. The shape-prediction accuracies were then analyzed by comparing the 
predicted values with the analytically-predicted values obtained from the finite-element analyses. The 
results showed that by using the piecewise nonlinear strain representations in the formulation of the 
Displacement Transfer Functions, the shape-prediction accuracies were greatly improved, especially for 
highly-tapered beam structures. 

NONUNIFORM AND UNIFORM BEAMS 

In the present technical paper, the term “beam” implies the imaginary embedded beam, which is 
defined as the depth-wise cross-section of the structure along the strain-sensing line. The “nonuniform 
beam” in the present technical paper is defined as the embedded beam with a varying depth factor, c(x) , 
a vertical distance from the neutral surface to the bottom surface of the nonuniform beam at axial 
location, x. For a nonuniform tubular beam, c(x) will be the outward radius. The “slightly nonuniform 
beam” in the present technical paper is defined as the embedded beam with slowly changing depth factor, 
c(x) . The “uniform beam” in the present technical paper is defined as the embedded beam with constant 
depth factor [that is, c(x) = c]. 

EMBEDDED BEAM CURVATURE EQUATIONS 

Figure 1 shows a deformed uniform embedded beam with constant depth factor, c. In a small beam 
segment bounded by dO , the small un-deformed curve length, AB{= ds ), is on the beam neutral surface, 

and the deformed curve length, A f B T {= AB[l + £(x)] = ds[\ + £(x) | } , lies on the beam bottom surface 

under bending strain, s{x). From the two similar slender sectors, O’AB and O t A t B t , the local radius of 
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curvature, of the beam elastic curve, can be related to the bottom surface bending strain, ;.'(x) , 

through the beam depth factor, c, as: 


O'A' 
O *A 


= 1 + 


A'B' 


R{ x) AB 


= 1 + £^JC) 


( 1 ) 


Applying the Lagrangian curvature equation, 1 / R(x ) = (d 2 y/dx 2 )/^l-(dy/dxf (refs. 12 and 

13), neglecting the term, (dy / dx) 2 (negligible axial displacements), equation (1) can be written in the 
form of a second-order differential equation as: 


d 2 y s{x) 
dx 2 c 


( 2 ) 


In equation (2), x is the undeformed axial coordinate, and y is the out-of-plane deflection. Because 
equation (2) is referred to the un-deformed x-coordinate (Lagrangian formulation), the term, d 2 yidx 2 , in 
equation (2) is not the simplified form of the classical curvature equation, 
\! R(x)=(d 2 y I dx 2 )!\Y + (dy I dx ) 2 ^ n , which is referred to the deformed x-coordinate (Eulerian 
formulation) (refs. 12, 13, and 14). 

Notice that curvature equation (2) is purely a geometrical relationship, containing no material or 
structural properties. In fact, the effect of the material or structural properties is felt by the surface 
bending strains, s{x)\ the value of s{x) is influenced by the material or structural properties. 

The curvature equation (2) could also be extended to the nonuniform embedded beam without losing 
accuracy if the embedded beam depth factor, c(x), changes slowly along the beam axis (ref. 15). 
Namely, 


cfy s(x) 
dx 2 c(x) 


( 3 ) 


For given beam geometries, the depth factors {c, c(x)} are known. If the functional form of the 
surface strain, s{x ), is specified, equation (2) or (3) can be integrated once to yield the beam slope 
equation, and the second integration will yield the beam deflection equation. To integrate the curvature 
equations (2) or (3), the piecewise integration method can be used. (refs. 1-5). 


DISCRETIZATION OF NONUNIFORM EMBEDDED BEAM 

The integrations of curvature equation (2) or (3) can be achieved by means of piecewise integrations 
through beam discretizations (refs. 1, 2, and 5). As shown in figure 2, the embedded beam of length, /, is 
first discretized evenly into n multiple small domains, x i X < x < x £ (/ = 1,2,3,...., n) , with equal domain 

length, A l(=l/n). The strain-sensing stations where the bending strains will be measured (or 
analytically calculated) are located on the beam bottom surface at the domain junctures (including at both 
ends of the beam), x = x 0 ,x 1 ,x 2 ,.„.,x n (fig. 2). Through discretization, both beam depth factor and 

surface bending strain, {c(x ) , £r(x)} , within each small domain, x it < x < X-, could be expressed with 
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simple linear functions. Such an approach enables the piecewise integrations of curvature equation (2) or 
(3) in closed forms for each small domain to yield beam slope and deflection equations in recursive 
forms. These recursive slope and deflection equations can be combined into a dual summation-form 
deflection equation called the Displacement Transfer Function (refs. 1-5) and expressed in terms of 
domain length, A/, beam depth factors, and surface bending strains, {c iy £ f } , (fig. 2). 

SUMMARY OF PREVIOUS DISPLACEMENT TRANSFER FUNCTIONS 

All the earlier Displacement Transfer Functions were formulated through piecewise integrations of 
equation (2) or (3). Four earlier Displacement Transfer Functions formulated for nonuniform, slightly 
nonuniform, and uniform straight embedded beams, are shown below for comparison of their 
mathematical functional forms (refs. 1, 2, and 5). 

1. Displacement Transfer Function for nonuniform beams (c t ^ c i Y ) (derivations in Appendix A of 
ref. 2): 


j-i 


S J-i~ S J 


£ J-i C J~ £ J C J - 1 


2(Cy_i Cj) (. c j-i C j ) 


6v]°g B j +{ Cji -Cj) 


J-l 


+(a/) 2 i; 

j = i 


(j-j) 


Contributions from deflection terms 


£jl £j + £j ' Cj £jC { ' ] 0 g e C J- 


(4) 




J - 1 


Contribution from slope terms 


+ y o + (i)Altaa0 o 


=0 for cantilever be ams 




2. Displacement Transfer Function for slightly nonuniform beams (c i / c i X — » 1) (derivations in 
Appendix C of ref. 2): 


y,=-r-Z 


J=1 


r 

1 

rr o 

(A/) 2 

i-1 

v i 

r 

c i-j ) 

rr c.'i "i 



T + - — — 
1 2 





(5) 


Contribution from deflection terms 

+ y 0 + A/tan 0 O 

' V ' 

=0 for c anti lever beams 


Contributions from slope terms 


0 = 1,2, 3, 

which was obtained from equation (4) by expanding the logarithmic terms in the neighborhood of 




3. Second-order Displacement Transfer Function for nonuniform beams ( c i ^ c i X ) [1 /c(x) in 
equation (3) expanded in series up to the second-order term] (ref. 5): 
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( 6 ) 


y= (Al)2 ± 1 1 
y ‘ 60 Ac . 

[ Cj r Cj VI 

27-9^ + 21 -L- | 

£ j - 1 + 

r r VI 

18 — 1 1 — — h 3l 

"1 

s i 

II 

vj 

1 

1 

\_ 

\ 


w. 1 

\_ 

/ 

\ 

> 


Contributions from deflection terms 


, (A/) 2 

y ( x A 

Cj f Cj VI 

Q — A J -u J l 

r 4- 

r r VI 

■N 

+ 12 

^ c 
>i j-i 

1 

15 

h 

1 

£ J 1 + 

1 

b 

b 

c 

V* 

1 

£ j 


Contributions from slope terms 


+ T Ko + (QA/tan<9 l) 
=0 for cantilever beams 


0 = 1 , 2 , 3 ,. ...,») 

4. Displacement Transfer Function for uniform beams ( c — c_ , = c ) (derivations in ref. 1; Appendix 
D of ref. 2; Appendix A of ref. 5): 


y , = 


(A /) 2 

6 c 


i-i 

(3/ - IK +6£(/ 

y=i 


+To +(/)A/tan6> 0 

V V ' 

=0 for cantilever beams 


0 = 1,2,3,...., «) (7) 


Equation (7) is the degenerated form of equations (5) and (6) for the uniform beam case, 
(C'=C'_ 1 =c). 

CHARACTERISTICS OF DISPLACEMENT TRANSFER FUNCTIONS 

Except for equation (4), equations (5) and (6) could degenerate into equation (7) for the uniform beam 
case (c L = c i X — c ) without incurring mathematical indeterminacy problems (refs. 1 and 2). Note that in 

each of the deflection equations (4)-(7), the deflection, at the current strain-sensing station, x v is 
expressed in terms of domain length, A/, inboard beam depth factors, (c 0 ,c 1? c 2 ,c 3 ,....,c f ), and the 
associated inboard surface strains, ( , s x , s 2 , s z e t ), including the values of ( c i , s t ) at the current 

strain-sensing station x 2 , where deflection, y t , is calculated. Because equations (4)— (7) contain no 

material or structural properties, when equations (4)— (7) are used in the deformed shape calculations of 
complex structures (for example, aircraft wings), there is no need to know the material or structural 
properties. In fact, the effect of the material or structural properties is felt by the surface strains, s t 

(*- 1 , 2 , 3,. ...,») . This is the key characteristic of the Displacement Transfer Functions (refs. 1-5). 


PREVIOUS PIECEWISE LINEAR REPRESENTATIONS 

The above four displacement transfer functions, equations (4)-(7), were formulated through 
discretization of the embedded beam into multiple small domains (fig. 2). The variation of the beam depth 
factor, c{x ) , and the surface bending strain, s{x ), within each small domain, x i X < jc < JC-, between the 

two adjacent strain-sensing stations, {x zl , x z } , was then expressed with linear functions as: 


c(x) = c i _ l +(c i -c i _ l ) 


x ~ x , , 

A/ 


(A , <*<*,) 


( 8 ) 
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e(x) = £ i _ l +(£ i -£ il ) 


CVi <x<x,) 


( 9 ) 


x — X 


A / 


Equation (8) can be a good piecewise linear representation of the slowly- varying depth factor, c(x) , 
of the nonuniform beams. The accuracy of equation (9) for piecewise linear representation of varying 
surface strain, s{x ), however, needs to be examined when the beam taper angle increases and causes 
strain distributions to be highly nonlinear. 


BEHAVIOR OF STRAIN CURVES FOR DIFFERENT TAPERED BEAMS 

This section is to examine the accuracy of piecewise linear strain representation equation (9) using 
tapered beams with different taper angles, and to show that equation (9) can lose accuracy in representing 
nonlinear strain distributions for the cases of highly-tapered beams. 


Tapered Cantilever Tubular Beams 

To study the accuracy degradation of the piecewise linear strain representations equation (9), a family 
of tapered cantilever tubular beams with different taper angles was chosen. The dimensions and taper 
angles, a , of the aluminum tapered cantilever tubular beams considered are listed in Table 1. 


Table 1. Dimensions of aluminum tapered cantilever tubular beams. 


l, in. 
(length) 

t , in. 

(wall thickness) 

c 0 , in. 
(root depth) 

in - 

(tip depth) 

c n /c 0 

(depth ratio) 

«r{=tan- 1 [(c 0 -c„)//]}, 
deg (taper angle) 

100.5 

0.02296 

4 

4 (uniform) 

4/4 (1.0) 

0.00 

100.5 

0.02296 

4 

3 

3/4 (0.75) 

0.57 

100.5 

0.02296 

4 

2 

2/4 (0.5) 

1.14 

100.5 

0.02296 

4 

1 

1/4 (0.25) 

1.71 

100.5 

0.02296 

4 

0.5 

0.5/4(0.125) 

1.99 

100.5 

0.02296 

4 

0.25 

0.25/4 (0.0625) 

2.14 


The imaginary embedded beam shown in figure 2 can also represent the embedded beam of a typical 
tubular beam, which is defined as the span-wise cross section through vertical diameters along the 
strain-sensing line. For each tapered tubular beam, the depth factor, c v of the embedded beam is the 

known local outer radius. Therefore, only one strain-sensing line on the bottom surface oriented in the 
span-wise direction is needed for shape predictions (fig. 2). 


Figure 3 shows the finite-element model of a typical tapered cantilever tubular beam (depth ratio, 
c n /c 0 = 2/4) generated from the SPAR (Structure Performance And Resizing) finite-element computer 

program (ref. 16). The size of the SPAR model (identical for all tapered beams) is also indicated in figure 
3. The tapered cantilever tubular beam is subjected to upward load of ,P=100 lb at the beam tip. 


Analytical Surface Strains 

In the present technical paper, the surface strains at the strain-sensing stations (fig. 2) were calculated 
analytically (not actually measured) using the SPAR program. Namely, finite elements were used to 
simulate actual strain gages. The surface bending strains, s t (/ = 0,l,2,3,....,w), at the i - th strain-sensing 
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stations (fig. 2) were obtained by converting SPAR-generated nodal or element axial stresses into axial 
strains, ^(bending strains), through stress-strain law. 

Also, the surface strain, , can be obtained from the span-wise nodal displacement differentials of 

the SPAR elements. It was found, however, that in the beam-tip regions of highly-bent tapered beams, the 
SPAR element strains became negative because the output nodal displacements were the projected 
displacements along the original beam axis, and not along the deformed beam axis to reflect true span- 
wise displacements. Therefore, this displacement method was not used. 

Shapes of Strain Curves 

The SPAR-generated strain curves were approximated with piecewise linear functions using domain 
densities of n = 8 and n — 16 . 

For n = 8 

Figure 4 shows the shapes of piecewise linear strain curves (solid lines with solid circular symbols) 
using domain density, n — 8, to approximate the SPAR-generated strain curves (dashed curves with open 
circular symbols) for different tapered cantilever tubular beams. Note that the SPAR-generated strain data 
are not exactly zero (theoretically zero) at the free beam tip (especially for the c n /c 0 = 0.5/4, 0.25/4 

cases); the reason being that the beam-tip finite element (simulated strain gage) has one end at the beam 
tip, but the other end is slightly off from the beam tip because of finite length. By reducing the element 
sizes, the beam-tip strain should approach zero. If an actual strain gage were installed at the beam tip, the 
strain output could also be nonzero because of finite length of the strain gage. 

For moderately-bent strain curves (for example, c n /c 0 = 3/4, 2/4, 1/4) (fig. 4), the piecewise linear 

strain representations could be reasonably good, but for highly-bent regions of the strain curves (for 
example, c n /c 0 = 0.5/4, 0.25/4), the piecewise linear strain representations could lose accuracy. Note 

from figure 4 in the highly-bent regions of the strain curves (for example, c n /c 0 = 0.5/4, 0.25/4), the 

typical piecewise linear points, calculated from equation (9), lie slightly off from the corresponding SPAR 
data points. 

Keep in mind that for the low depth ratios, (for example, c n /c 0 = 0.5/4, 0.25/4), each beam will bend 

like a fishing pole in the outboard regions. These unrealistic very low depth ratios were included in the 
current analysis solely for theoretical interest to examine the shape-prediction accuracies of the 
Displacement Transfer Functions. For real aircraft wings, the depth ratios are much higher. For example, 
the Ikhana wing has a depth ratio of c n /c 0 = 0.4444, for which the strain curve is bow-shaped (ref. 3). 

For n = 16 

To improve piecewise linear strain representations using equation (9), the domain density was 
increased from « = 8 to n — 16 to cause the piecewise linear strain curves to approach the corresponding 
SPAR-generated smooth strain curves. 

Figure 5 shows that by using domain density n — 16, the piecewise linear strain representations 
greatly improved except for the sharply-bent regions. The peak strain point at strain-sensing station 
/ = 15 , which was missed at « = 8 for the cj c 0 = 0.25/4 case, is now included. In the sharply-bent 
region of the strain curve of the c n /c 0 = 0.25/4 case, the typical piecewise linear point calculated from 
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equation (9) lies below the corresponding SPAR data point; however, further increasing of the domain 
density, n , could cause the piecewise linear strain curves to approach the corresponding SPAR-generated 
smooth strain curves. 

An alternative approach is to use piecewise nonlinear strain representations (described below) to 
approximate the nonlinear strain curves and thereby improve the shape-prediction accuracies. 

FORMULATIONS OF IMPROVED DISPLACEMENT TRANSFER 

FUNCTIONS 

In the formulation of the Improved Displacement Transfer Functions, the piecewise linear assumption 
according to equation (9) for strain distribution was replaced with the piecewise nonlinear strain 
representations described below. Equation (8) was used without modification, however, because 
equation (8) is a good piecewise linear representation of the slowly- varying depth factor, c(x) , of the 
nonuniform beams. 


Piecewise Nonlinear Strain Representations 

In the formulation of the Improved Displacement Transfer Functions, the distribution of surface 
bending strain, s{x ), within the dual domain, x iX <x< X l+l , was described with piecewise nonlinear 
(second-order) functions in terms of three strain values, { s iX , , £ i+l } , respectively at three adjacent 

strain-sensing stations, { x zl , X l , x +1 }, (fig. 6). The resulting piecewise nonlinear strain equation has the 
following form: 


4 *) = 1 


3s . , - 4s. + s M/ . s, -2s + s i+ . . . 2 


2A/ 


2(A If 

(j = 0,l,2,3,....,w) 


( 10 ) 


Equation (10) is the piecewise nonlinear strain equation to replace the piecewise linear strain equation 
(9) for formulating the Improved Displacement Transfer Functions. 


Strain Extrapolation 

Note that equation (10) is for the strain-sensing stations x i =x 0 ,x 1 ,x 2 ,x 3 ,....,x n . At the beam tip 
strain-sensing station x t = x n , however, an extra strain point, s n+l , is needed for three-points strain curve 
approximation. A simple way to estimate s n+l is through extrapolation using the three known strains, 
respectively at the three adjacent strain-sensing stations, {x n _ 2 ,x n 1? x n }, (fig. 7). The 
three-point extrapolation equation established for calculating s n+1 has the following form: 


'n\ 1 


= 


n- 2 


- 3r. . + 3r. 

n\ n 


( 11 ) 


Equation (11) was obtained from equation ( 1 0) by setting i = n — 1 and x — x j = x w|| — x n 2 = 3 A/ . 
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Piecewise Nonlinear Strain Curves 

Figure 8 compares the piecewise nonlinear strain curves (solid curves with solid circular symbols) 
based on equation (10) using the domain density n = 16, and the corresponding SPAR-generated strain 
curves (dashed curves with open circular symbols) for different tapered cantilever tubular beams. In the 
sharply-bent segment of the strain curve of the c n /c Q = 0.25/4 case (fig. 8), the typical piecewise 

nonlinear point, calculated from equation (10), has now moved above the corresponding SPAR data point. 
Except for the c n /c 0 = 0.25/4 case, the piecewise nonlinear strain representations based on equation (10) 

could provide very good approximations of the bow-shaped SPAR strain curves. As will be seen, by 
changing from piecewise linear to piecewise nonlinear strain representations (see figs. 5 and 8), shape- 
prediction accuracies could be greatly improved. 

UNIFORM BEAMS 

The formulation of the Improved Displacement Transfer Functions will start from the mathematically 
simpler case of uniform beams, and progress to the mathematically more complex case of nonuniform 
beams. As mentioned above, the uniform beam considered in the present technical paper is the embedded 
beam with constant depth. In the formulations of the Displacement Transfer Functions, dimensions of 
beam width are not needed. Using the piecewise linear depth equation (8), and the piecewise nonlinear 
strain equation (10), the Improved Displacement Transfer Functions were formulated first for the uniform 
beams. The slope and deflection equations are to be obtained by piecewise integrations of the curvature 
equation (2) for the uniform beams. 


Improved Slope Equation 

The slope, tan#(x)D of the uniform beam at axial location xQvithin the small domain, 
x i X < x < X i? bounded by the two adjacent strain-sensing stations, l , x i } , can be obtained by 

integrating the uniform-beam curvature equation (2) once, and enforcing the continuity of slope at the 
adjacent inboard strain-sensing station, x i X . Namely, 


tan*9(x) = f —^rdx + taxi0 i X = f *^-^-dx + tan#^ ; (jc < x < x ) 
J vi dx v — v — ' J *i i c 


s{x) 


Slope increment 


Slope at x / _ ] 


C 

Eq. (2) 


( 12 ) 


in which tan0 , is the slope at x ; , . Substituting equation (10) into equation (12) and carrying out the 
integration (ref. 17), one obtains the slope equation for the small domain, jc x <x<x t (see Appendix A): 


tan 0(x) = 


£ n( x - x n)- 


3c t , -4c, + c, 


2+1 


4A/ 


(x-x,) 2 + 


'i - 1 


2c, + £, 


2+1 


6(A/) 2 


(x-x | ) 3 


+ tan6>_, (13) 


(Xj_i <*<*,) 
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At the strain-sensing station, x — x , we have x i —x j _ l =Al, and equation (13) gives the slope, 
tan 0\= tan 0{x i )] ,at the strain-sensing station, x v as: 

tan = 7 ”X 5 ^-i + “ £ m ) + ^ 3-i ; (i = l2,3,....,n) (14) 

12c 

Equation (14) is the Improved slope equation for the uniform beams in recursive form. Applying the 
recursive relationship, equation (14) can be written as: 

A/ 1 

fan0i = — Z( 5 ^-i +8 ^-^ + t)+tan0 o ; (/ = 1,2,3,...., n) (15) 

lzc j—\ 

Equation (15) is the Improved slope equation for the uniform beams in summation form. 

Improved Deflection Equation 

The deflection, y(x ), of the uniform cantilever beam within the small domain, x iY < x < can be 

obtained by integrating the slope equation (13) and enforcing the continuity of deflection at the inboard 
adjacent strain-sensing station, jc- _ 1? as: 


X*) = £ i Jan 0(x)dx +y^ ; (x M < x < x, ) 

eq. (13) 


(16) 


in which y i { is the deflection at the inboard strain-sensing station, jc z1 . 

In view of equation (13), equation (16) can be integrated (ref. 17) to yield the following deflection 
equation for the small domain, l < jc < x n as (see Appendix A): 


>«=- 

c 


^-(x-x,-,) 2 - ^ 1 (X - + *»„*? *^ (x-x,,,) 4 


24(A/) 


(17) 


+(*-* i _i)tan0 M +y._ l 


(x, i <x<x,) 


At the strain-sensing station, x = x n we have x t — x i V = A/, and equation (17) gives the deflection, 
, at the strain-sensing station, x ; , as: 


y, 


(A If 

24c 


(7^, + 6^, - +y i _ 1 + M tan 0^ ■ (/ = 1,2,3,....,/?) 


(18) 


Equation (18) is the Improved deflection equation for the uniform beams in recursive form. 
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Improved Displacement Transfer Function 

In view of equation (15), and applying recursive relationships, equation (18) can be written in dual 
summation form as (see Appendix B): 


y, = ^ 7 Z ( 7 ^1 + 6 - £ J+l )+y 0 + X (* - ! + 8 ^ -^ + i)+0)Vtan<9 0 (19) 


Contributions of deflection tenns 


Contributions of slope terms 


(i-1.2,3,— -,«) 


Equation (19) is called the Improved Displacement Transfer Function for the uniform beams because 
equation (19) transforms the surface strains into deflections. Note that the first summation is contributed 
from the deflection terms, and the second summation is contributed from the slope terms. Note also that a 
factor (i—j) appeared in the second summation term due to cumulative summations of identical terms 
(see Appendix B). Note that the second summation ends at j = i — \ (that is, one term shorter than the 
first summation) (see Appendix B). 


NONUNIFORM BEAMS 

As defined earlier, the nonuniform beam has a variable depth factor, c{x) . Using the piecewise linear 
depth equation (8), and the piecewise nonlinear strain equation (10), the Improved Displacement Transfer 
Functions were formulated for the nonuniform beams. The slope and deflection equations are obtained by 
piecewise integrations of the curvature equation (3) for the nonuniform beams. 

Improved Slope Equation 

The slope, tan #(:*:)□ of the nonuniform beam at axial location x within the small domain, 
x t l < x < x i? bounded by the two adjacent strain-sensing stations, {x zl , x £ } , can be obtained by 

integrating the nonuniform-beam curvature equation (3) once, and enforcing the continuity of slope at the 
adjacent inboard strain-sensing station, x t l . Namely, 

“ n = e * + ,an 5-i = c ian ^ ■ ; a '"- isx ' ) <20) 

v v , Slope at * v t 

Slope increment Eq. (3) 

in which tan#. 1? is the slope at x t l . Substituting equations (8) and (10) into equation (20), and carrying 
out integration, one obtains the slope equation for the small domain, x i X < x < x r as (see Appendix C): 
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A/ 


tan6»(x) = — -[(2^-c^Xc^ -2c I _ l f J ) + c J c i _ 1 ^ +I ]x 


2(^-1) 


1 

2(^-c l _ 1 ) : 


X 1 log. 




x-x. 


A/ 


-log^,, 


[(2c - -1c, , K, + 2 (2c - Cj_j )c - a;, ] (x - x, , ) 


( 21 ) 


+ g , 1 -2 gi + gm (x-x il ) +ta ^ 
4 (c M 


(x ;1 <x<x,) 


At the strain-sensing station, x = x p we have x, — x (1 =A/, and equation (21) gives the slope, 
tan$ ; L= tant)(x)| , at the strain-sensing station, x t , as: 


tan 0 = 


A/ 


' 2 ( c ,- c ,_ 1 ) 3 

A/ 

4( c , _ Ai) 2 


[(2c, -c^Xc s , -2c, l ^) + cc_,i- +l ]log B -L. 

C i 1 

[(5c, -3c, -2(3c, - c, t )*, + (c, + c M )ff i+1 ] + tan <9 , 

0 = 1 , 2 , 3 ,. 


( 22 ) 


Equation (22) is the Improved slope equation for the nonuniform beams in recursive form. Applying the 
recursive relationship, equation (22) can be rewritten as: 


tan6 '< = vE 7 1 — a[( 2c , - c )(c c - 2c e.) + c,cc]log e Cj - 

2 j^[( c j~ c j-i) J c j-i 

~ 2(c -c ) 2 [ (5c J - 3 c j -.H-i“2(3c j -Cm>v + (c y + c y _ I )c /+I ]| + tan6> 0 

0-1,2,3,....,k) 


(23) 


Equation (23) is the Improved slope equation for the nonuniform beams in summation form. 

Improved Deflection Equation 

The deflection, y(x ), of the nonuniform cantilever beam within the small domain, x i X < x < jc-, can 

be obtained by integrating the slope equation (21), and enforcing the continuity of deflection at the 
inboard adjacent strain-sensing station, x i X , as: 


A x ) = J +y f _ , ; (x, t < x < x, ) 

eq. ( 21 ) 


(24) 
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in which y i X is the deflection at the inboard strain-sensing station, x t l . 

In view of equation (21), equation (24) can be integrated (ref. 17) to yield the following deflection 
equation for the nonuniform beam for the small domain, x iX < x < x r as (see Appendix C): 


y(x)= 


Al 


2{c- Ci J 
xlogjc, +(c,-c,_,) 


[( 2c , - C i-i Xc,*,- 1 - 2c i-i s i ) + C i c i-S ; M ] 

( c I“ c I-l) 


Al 


( c I + (c | .-c i ._ 1 ) 


x-x t 

Al 


Al 


Al 


C i-l lo Se C /-l )-(*- Vl) lo g* C ,1 


(25) 


+ -r— 1 — ^[( 2c > ~ 3c mKi + 2 ( 2c i - - c i £ i+ i]( x ~ x i- 1) 2 

4(c, -c^) 

^,- 2 ^+^ ( x - x ,_ t ) 3 

+ —^7 7 . , + y^-i + (* - Vi ) to 11 

12( C| -c M ) A/ 


(x (1 <x<x,) 


At the strain-sensing station, x = x n we have jc, — a, , = A l , and equation (25) gives the deflection, 
yj= y(-t )| at the strain-sensing station, x r as: 


y,= 


(Alf 


(V ) 2 

1 2(c- Ci J 


[( 2c i C i-l)( C i £ i-l 2c /-l £ j) + C l C i-l £ i + i] 




(26) 


[(8c, -5c ,)c- , - 2(5c -2c, ,)&■ +(2c,+c, ,)£.', ,]+_v, , +A/tanO, , 


0 = 1,2,3,....,«) 

Equation (26) is the Improved deflection equation for the nonuniform beam in recursive form. 


Improved Displacement Transfer Function 

In view of equation (22), and applying recursive relationships, equation (26) can be written in dual 
summation form as (see Appendix D): 
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<7 log . — ) 


X=(A/) ! £ 

J=1 


2 ( 6 --c- ,y t (2c -i - ^ c j-i s j ) + c j c j-\ £ j+\ ] 

- i2( c ! C )2 [(*7 - 5c m) £ m - 2 ( 5c ; - 2c y-i) g j + ( 2c y +S-i)v] 


•+y 0 


Contribution from deflection terms 


J = 1 


— — [(2 Cj -c, , X^^-, - ^ c j-i £ j ) + c j c j-i £ j+ i ] 1°g B 
z V c y c y-iJ c j-i 

4 , — ~^[< 5 S -3c 1 )£ l -2(3c -cjs +(c +C ,)^ +1 ] 

^ c y c j-ri 


>+(/)A/tan6? 0 


Contributions from slope terms 


(z = U,3,....,«) (27) 


Equation (27) is called the Improved Displacement Transfer Function for the nonuniform beam. 
Notice that the first summation is contributed from the deflection terms, and the second summation is 
contributed from the slope terms. 


LOG-EXPANDED CASE 

Note that the slope equations (22) and (23) and the deflection equations (26) and (27) contain 
logarithmic terms, log e {cJ c i X ), in the numerators, and (c i — c i V ) terms in the denominators, causing 

mathematically indeterminacy (0/0) when the beam depth becomes uniform (that is, c % — c i X — c). This 

mathematical indeterminacy problem can be eliminated if the logarithmic terms are expanded in terms of 
(Pi-*- 1) series to cancel the ( c i — c i X ) terms in the denominators (refs. 1, 2). 


Log-Expanded Slope Equation 

For (c E / c^)— » 1, the logarithm term, log^c, / c^), in equations (23) and (24) can be expanded in 
series form up to the third-order term as (ref. 17): 


, c, (c, t ) i ( c, y i tc, y 
^•yyUr'n lyy-'J 3lyy _1 J -• 

= +(c, -c,.,)<2c, - 5c,.,)-....] 


(28) 


Substitution of equation (28) into equations (22) and (23), and after grouping terms, yields the 
Log-expanded slope equations given below (see Appendix E): 


15 


1 . In recursive form: 


tan 3 = {[5c? , + 4{c, - c M ) 2 ]^., - 8c, .^(c, - 2c,_ 1 >.- + c,. t (2c, - 3c, ,)c w } + tan6>, (29) 

(* = 1.2,3,....,/?) 

2. In summation form: 


ton ^ = Z ^ {[ 5c > 1 + 4 ( c j - C J - 1 ) 2 1^-1 - 8c m ( c j - 2c M H- +c ,-i (2Cj - 3c ,)r„ } 

12 j=i c j i 


(30) 


+ tan ft. 


0 = 1,2,3,....,/?) 


Log-Expanded Deflection Equation 

For developing the Log-expanded deflection equations for (c, / c, , ) — > 1 . the logarithm term, 
1 og e (c, / c, ,). in equations (26) and (27) must be expanded in series form up to the fourth-order term as 
(ref. 17): 


1 C i ( C i 3 1 ( c t ^ if c t ^ if c, ^1 

= ^H 12 4, -6c*,(c, -c,,)+4c,.,(c, -e,_,f-3((c, -c,J +....] 


( 31 ) 


Substitutions of equation (31) into equations (26) and (27), and after grouping terms, yields the 
resulting Log-expanded deflection equations given as (see Appendix E): 


y l 


(A If 

24c, 4 , 


[7c? j + (c, - Cj.jXSC/Cj-.j - 3c?>- 3(c? -cl)] £i l 
+2c, , [(3c? + (c, - c, , )(3c, - 8c, t )]c - c, , [c? + 2(c f - c, , f ] c, ( 

+_p ; | + A/ tan ft, , 

(/ = 1,2,3,....,/?) 


(32) 


Equation (32) is the Log-expanded deflection equation for the nonuniform beam in recursive form. 
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Log-Expanded Displacement Transfer Function 

In view of equation (29), and applying recursive relationships, equation (32) can be written in dual 
summation form as (see Appendix F): 


(A Iff 1 , 
24 UcU 

V 

K- 

+*>.i 

- 3cj)- 3(, 


h- 

y J+ n Cj - Cj j] 


v 

+ 

* 

(A/) 2 y (/-y) \\_^ C j 
12 £ c)_, [ + c,_ 

Contributions from deflection terms 

_! + 4 <c - c _ t ) 2 ~\s , - 8c (c - 2c ! )c 

J J J J J J J J U(/)A/tan<9 0 

i (2 C j ~ 3 Cj x )c J+1 

(33) 


Contributions from slope terms 

(/ = 12 , 3,. ...,«) 


Equation (33) is called the Log-expanded Displacement Transfer Function for nonuniform beams. 


Log-Expanded Slope Equation for Uniform Beams 

For uniform beams, we have c i — c iX — c , and the above Log-expanded slope equations (29) and (30) 
degenerate respectively into the following simplified forms: 

1 . In recursive form: 


tan 0 i = — (5c, , + 8^ - , ) + tan 0 t , ; (/ = 1,2,3,....,/?) (34) 

12c 

2. In summation form: 

tan 0, = ^ L(5c, , + 8c, - c, , ,) + tan 0 O ; (/ = 1,2,3,....,/?) (35) 

y=l 

Equations (34) and (35) agree respectively with equations (14) and (15) formulated for the uniform 
beams, and thus confirm the derivation accuracy of the Log-expanded slope equations (29) and (30). 

Log-Expanded Deflection Equation for Uniform Beams 

For uniform beams, we have C- = c iX — c , and the above Log-expanded deflection equations (32) and 
(33) degenerate respectively into the following simplified forms: 

1 . In recursive form: 

y i = l^(7 £ i-i+ 6£ i- £ M )+yi-i+Mtato0 i ^ ; (/ = 1,2,3,....,/?) (36) 
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2. In summation form: 


(Klf 1 (\l\ l 11 

y< = +6s j ~ £ j+i) + yo +g£ j -£ j«)HQM tan6> 0 (3y) 


Contributions from deflection terms 


Contributions from slope terms 

0 = 1 , 2 , 3 , 


Note that equations (36) and (37) agree respectively with equations (18) and (19) formulated for the 
uniform beams, and thus confirm the derivation accuracy of the Log-expanded deflection equations (32) 
and (33). 


DEPTH-EXPANDED CASE 

This section describes the formulation of the Depth-expanded slope and deflection equations for 
nonuniform beams. The purpose of this formulation is to avoid the logarithmic terms after carrying out 
direct integration of curvature equation (3). The depth factor of the nonuniform beam is assumed to 
change slowly along the beam axis. 


Depth-Expanded Slope Equation 

The slope, tan#(x)D of the nonuniform beam at axial location, x, within the small domain, 
x ir <x<x r between the two adjacent strain-sensing stations, {x i{ , x z } , can be obtained by carrying 
out the integration of the nonuniform beam curvature equation (3). Namely, 

t an^AT ) = J'^A+toi^ : ,=£ i ^|<fc + tan^ 1 ; (*,_,<*<*,) (38) 

v v ' Slope at x, j i — v — j 

Slope increment pq (3) 

in which tan 0 i X is the slope at x z l . 


An alternative way to avoid the logarithmic terms after direct integration of equation (38) is to expand 
the term, 1 / c(x) in the integrand of equation (38) in binomial series form because c{x) is a 
slowly- varying function of x (ref. 17). Namely, 


1 

c(jc) 



(39) 


Note that in equation (39), the variable term (jc — x L l ) in the denominator was moved to the 
numerator. This approach avoids the appearance of logarithmic term, log^cv / c iX ), after direct 
integration of equation (38). Substituting equations (10) and (39) into equation (38) and carrying out 
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integration, one obtains the Depth-expanded slope equation for the small domain, x iX <x<x t , as (see 
Appendix G): 


tan 0(x) = <j £ iX (x — x lX ) - 


3^i- 4 ^ +%i 


+ £ 


( C ^ 

itr 1 ) 


2A1 


+ 


c . -2c+c . 3c f . -4c + c j+1 f c, ( c t ,^l 

' 1 1 '"+ — '--lj + c,_ t |^— -lj 


(x-x, ,) 3 


3(A0 2 


(40) 


G 1 2^ + £T /+1 f c t 

d , 3c,_ 1 -4c ( +c 1+1 f c, 

.T1 

(x-x, ,) 4 

2 U_i 

J 2 U„ 

j 

4(A/) 3 


►+tan# 


(X, L < X < X ; ) 

At the strain-sensing station, x = x ; , we have x j —x j _ l =Al, and equation (40) gives the slope, 
lan#J= tan#(x f )] at the strain-sensing station, x t , as: 


tan 0, 


Al 


' 24 c. 


( c \ 

2(5^1 + - £) + i) - ^ - lj (3^-i + 10c, - c, +1 ) 


(41) 


( V 

-|^7 — lj (£,! -12^+3^,) + tan 

0 = 1,2,3,....,«) 

Equation (41) is the Depth-expanded slope equation for the nonuniform beams in recursive form. 

For the uniform beams, we have c i — c iX = c , and equation (41) degenerates into the following 
simplified forms: 


tan# 


M 


(5 c,_, + 8c, - c i+1 ) + tan 0 l , ; (i = 1, 2, 3, ... «) 


' 12c 

which agrees with equation (14) formulated for the uniform beams. 

Applying the recursive relationship, equation (41) can be rewritten as: 


(42) 


tan0. = — £ — 

24^ c 

^ j=i c j -l 


r c ^ 

2(5^_j + 8c - c j+1 ) - - lj (3c, + 10c - c +1 ) 

r c v 

t^T _1 J ( v.- ,2 ^ + 3^ + . ) 


(43) 


+ tan#. 


0 = 1, 2, 3,. 


19 


Equation (43) is the Depth-expanded slope equation for the nonuniform beams in summation form, 
and which can be degenerated into equation (15) for uniform beams = c i X — c). 

Depth-Expanded Deflection Equation 

The deflection, y(x ), of the nonuniform cantilever beam within the small domain, x i X < x < JC-, can 

be obtained by integrating the slope equation (40), and enforcing the continuity of deflection at the 
inboard adjacent strain-sensing station, x iX , as: 


yO) = { x tzn00c)dx+ y^ ; (x M < x < x,) 


Eq. (40) 


Deflection 
1 


(44) 


in which y i X is the deflection at the inboard strain-sensing station, x i X . In light of equation (40), 

equation (44) can be integrated (ref 17) to yield the following deflection equation for the small domain, 
x i X <x<x r (see Appendix G): 


, . i (x-x^y 

}<x) = — — — 

c,-, I 2 


3s t , -4e t +£ M ( C, ,1 

— - ' — — -i-^J — — 1 


itr 1 ) 


(x-x 4 1 ) 3 


6A/ 


£: . - 2 s + e 3s , . - 4s, + s ( c, ( c, . V 

11 1 ' +1 + — — L lj 




12 (A If ( 45 ) 


' i 1 i 1 1 1 i 


-H+ : 


Ur'J 

+y, ^(x-x^tan^ ! 

At the strain-sensing station, x = x,, we have x, — x ;l =Al, and equation (45) gives the deflection, 
V ; |=v(.v)| at the strain-sensing station, x,. as: 


K-iT" 


U , J _ 

20(A/) 3 


(x,i<x 


y-,= 


(Aiy 


24c, 


2 ( c } 

(? £ i X + 6 s, lj (4 + 7 s t - s M ) 

+ 5l^7 _1 J ( e « +12e i- 3ff «) 


(46) 


+y, , + A/ tan 0 i X 

0 -1,2,3,...., «) 

Equation (46) is the Depth-expanded deflection equation for the nonuniform beams in recursive form. 
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For the uniform beams, we havecr —c^ — c , and equation (46) degenerates into the following 
simplified form: 


* = -,+6c, ~e m )+y, .,+A/tanf,., ; (i = l,2,3,....,») (47) 

Equation (47) agrees with the deflection equation (18) for the uniform beams. 

Depth-Expanded Displacement Transfer Function 

In view of equation (41), and applying recursive relationships, equation (46) can be written in dual 
summation form as (see Appendix H): 


24 £< 7 -, 


2(c, ^ 

(7 s j-i + 6s j ~ £ j + 1)- 5 ^-!j (4^ i + 7s j- £ j+ 1) 

ir Cy v 


+ v„ 


24 & <v_, 


Contribations from deflection terms 

(c. } 


2(5£j_i +8£ J -£ J+1 )-\-*—lj(3£ J _ t +10 Sj -S J+1 ) 
( c V 

( £ j-i _ * 2 £+ 3 ^ +1 ) 


+(/)A/tan^ 0 


(48) 


Contributions from slope terms 


(z = l,2,3,....,«) 


Equation (48) is the Depth-expanded Displacement Transfer Function for the nonuniform beams, and 
which can degenerate into equation (19) for the uniform beams = c i X — c). Again, the first summation 

in equation (48) is the contributions from the deflection terms, and the second summation is the 
contributions from the slope terms. 


NEW STRUCTURE SHAPE-SENSING TECHNOLOGY 

The Displacement Transfer Functions with the accompanying strain-sensing system thus formed a 
new technology for structure deformed shape predictions. This new technology is called Method for 
Real-Time Structure Shape-Sensing (U.S. Patent Number 7,520,176, ref. 9), which could be very useful 
for in-flight deformed shape monitoring of flexible wings and tails, such as those often employed on 
unmanned flight vehicles by the ground-based pilot for maintaining safe flights. 

A typical flow chart of the new structure shape-sensing technology using the two-line strain-sensing 
system is shown in Appendix I. In the flow chart, only the Depth-expanded Displacement Transfer 
Function (48) was used as the example for each strain-sensing line. Other Displacement Transfer 
Functions may also be selected, depending on the types of structures (uniform, nonuniform, slightly 
nonuniform, et cetera) under consideration. 
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The two-line strain-sensing system can be used for the shape-sensing of the tapered wing box for 
which the neutral surface is located at the wing box half-depth (fig. 1-1), If, however, the neutral surface is 
not located at exactly half the depth of a structure, such as aircraft wings, a second set of the two-line 
sensing system must be installed on the upper surface, right on top of the lower sensing system. Using 
pairs of lower and upper strains, one can determine the actual location of the neutral surface for finding 
the values of the depth factors, c i (/ = 0 , 1 , 2 , 3 needed for the calculations of deflections (see the 
Ikhana wing shape predictions, ref. 3). 

It is important to mention that when using the two-line strain-sensing system for the wing structures, 
the Displacement Transfer Functions can not only calculate the deflections but also the cross-sectional 
rotations. Thus, the need to use surface distortion strains can be eliminated. 

Appendix I also describes the application of the two-line strain-sensing system to a free-free beam 
structure such as an airborne aircraft fuselage (free-free nonuniform tubular beam structure) (fig. 1-2) for 
sensing the horizontal and vertical bending strains for inputs to the Displacement Transfer Functions for 
the calculations of in-flight deformed shapes in horizontal and vertical planes (ref. 2). 

SHAPE-PREDICTION ACCURACIES 

As mentioned above, the surface strains were analytically calculated from the SPAR program, and 
were not directly measured. In the shape-prediction accuracy studies, the deflections calculated from the 
SPAR program were used as reference yardsticks (0% error) to evaluate the accuracies of the 
corresponding deflections calculated from the new Displacement Transfer Functions. The 
shape-prediction accuracies of the newly-formulated Displacement Transfer Functions were also 
compared with those of the earlier Displacement Transfer Functions. 

Tapered Cantilever Tubular Beams 

The sample structures chosen in the shape-prediction accuracy studies were again a family of 
aluminum tapered cantilever tubular beams having the dimensions listed in table 1. The same loading 
condition of upward load, P = 100 lb, was applied at each beam tip. 

Comparisons of Predicted Deflections 

Tables 2(a)-2(f) compare the accuracies of deflections calculated from different Displacement 
Transfer Functions based on piecewise linear and piecewise nonlinear strain representations for the six 
tapered cantilever tubular beams listed in table 1. Note from tables 2(a)-2(f) that only two typical 
Displacement Transfer Functions [the nonlinear Displacement Transfer Function (4); and the second- 
order Displacement Transfer Function (6)], based on piecewise linear strain representations were added to 
the accuracy comparisons. In tables 2(a)-2(f), the shaped prediction errors (in percent) are indicated in 
parentheses. 
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Table 2(a). Comparisons of deflections (in inches) calculated from SPAR and from different 
Displacement Transfer Functions; uniform cantilever tubular beam; c n /c 0 = 4/4; n = 16 . 


Theories 

yo 

Y2 

Y4 

Y6 

ys 

Yio 

y 12 

yi4 

yi6 

(% error) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01787 

(0.0000) 

0.06340 

(0.0000) 

0.13368 

(0.0000) 

0.22450 

(0.0000) 

0.33181 

(0.0000) 

0.45174 

(0.0000) 

0.57933 

(0.0000%) 

0.71162 

(0.0000) 


Piecewise linear strain representations 


Nonuniform, 
eq. (4)* 










Second- 

order, 

eq. (6) 

0.00000 

(0.0000) 

0.01576 

(0.2965) 

0.06034 

(0.4300) 

0.12962 

(0.5705) 

0.21948 

(0.7054) 

0.32580 

(0.8446) 

0.44448 

(0.9752) 

0.57139 

(1.1158) 

0.70243 

(1.2914) 


Piecewise nonlinear strain representations 


Improved, 
eq. (27)* 










Log- 

expanded, 
eq. (33) 

0.00000 

(0.0000) 

0.01577 

(0.2951) 

0.06036 

(0.4272) 

0.12965 

(0.5663) 

0.21951 

(0.7012) 

0.32585 

(0.8375) 

0.44454 

(1.0118) 

0.57146 

(1.1059) 

0.70249 

(1.2830) 

Depth- 
expanded, 
eq. (48) 

0.00000 

(0.0000) 

0.01577 

(0.2951) 

0.06036 

(0.4272) 

0.12965 

(0.5663) 

0.21951 

(0.7012) 

0.32585 

(0.8375) 

0.44454 

(1.0118) 

0.57146 

(1.1059) 

0.70249 

(1.2830) 


* No solutions for uniform beams because of mathematical indeterminacy. 


Note from Table 2(a) that both equations (4) and (27) produced no solutions for uniform beams 
because of mathematical indeterminacy problems. 

Table 2(b). Comparisons of deflections (in inches) calculated from SPAR and from different 
Displacement Transfer Functions; tapered cantilever tubular beam; c n /c 0 = 3/4; n — 16 . 


Theories 

Yo 

Y2 

Y4 

Y6 

ys 

Yio 

yi2 

yi4 

yi6 

(% error) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01823 

(0.0000) 

0.06697 

(0.0000) 

0.14572 

(0.0000) 

0.25204 

(0.0000) 

0.38294 

(0.0000) 

0.53428 

(0.0000) 

0.70092 

(0.0000) 

0.87634 

(0.0000) 


Piecewise linear strain representations 


Nonuniform, 
eq. (4) 

0.00000 

(0.0000) 

0.01627 

(0.2237) 

0.06417 

(0.3195) 

0.14199 

(0.4256) 

0.24738 

(0.5310) 

0.37726 

(0.6482) 

0.52758 

(0.7645) 

0.69312 

(0.8901) 

0.86711 

(1.0532) 

Second- 
order, 
eq. (6) 

0.00000 

(0.0000) 

0.01627 

(0.2237) 

0.06417 

(0.3195) 

0.14199 

(0.4256) 

0.24738 

(0.5310) 

0.37726 

(0.6482) 

0.52759 

(0.7634) 

0.69313 

(0.8888) 

0.86712 

(1.0521) 


Piecewise nonlinear strain representations 


Improved, 
eq. (27) 

0.00000 

(0.0000) 

0.01627 

(0.2237) 

0.06419 

(0.3172) 

0.14204 

(0.4199) 

0.24748 

(0.5203) 

0.37741 

(0.6310) 

0.52781 

(0.7383) 

0.69345 

(0.8524) 

0.86754 

(1.0042) 

Log- 

expanded, 
eq. (33)* 

0.00000 

(0.0000) 

0.01627 

(0.2237) 

0.06419 

(0.3172) 

0.14205 

(0.4188) 

0.24748 

(0.5203) 

0.37742 

(0.6299) 

0.52782 

(0.7372) 

0.69346 

(0.8513) 

0.86755 

( 1 . 0030 )* 

Depth- 
expanded, 
eq. (48) 

0.00000 

(0.0000) 

0.01627 

(0.2237) 

0.06419 

(0.3172) 

0.14204 

(0.4199) 

0.24748 

(0.5203) 

0.37741 

(0.6310) 

0.52781 

(0.7383) 

0.69345 

(0.8524) 

0.86754 

(1.0042) 


* Most accurate prediction (closest to SPAR). 
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Table 2(c). Comparisons of deflections (in inches) calculated from SPAR and from different displacement 
transfer functions; tapered cantilever tubular beam; c n /c Q = 2/4; n = 16 . 


Theories 

yo 

Y2 

Y4 

Y6 

ys 

YlO 

y 12 

Yl4 

yi6 

(% error) 

SPAR 

0.00000 

0.01864 

0.07113 

0.16057 

0.28831 

0.45510 

0.65965 

0.89759 

1.15801 

(reference) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

Piecewise linear strain representations 

Nonuniform, 

0.00000 

0.01682 

0.06859 

0.15719 

0.28407 

0.44983 

0.65325 

0.88981 

1.14828 

eq.(4) 

(0.0000) 

(0.1572) 

(0.2193) 

(0.2919) 

(0.3661) 

(0.4551) 

(0.5527) 

(0.6718) 

(0.8402) 

Second- 

order, 

eq. (6) 

0.00000 

0.01682 

0.06859 

0.15719 

0.28407 

0.44983 

0.65324 

0.88980 

1.14826 

(0.0000) 

(0.1572) 

(0.2193) 

(0.2919) 

(0.3661) 

(0.4551) 

(0.5535) 

(0.6727) 

(0.8420) 

Piecewise nonlinear strain representations 

Improved, 

0.00000 

0.01682 

0.06861 

0.15723 

0.28417 

0.45002 

0.65359 

0.89042 

1.14930 

eq. (27) 

(0.0000) 

(0.1572) 

(0.2176) 

(0.2884) 

(0.3575) 

(0.4387) 

(0.5233) 

(0.6192) 

(0.7522) 

Log- 

expanded, 
eq. (33)* 

0.00000 

0.01682 

0.06861 

0.15723 

0.28417 

0.45002 

0.65359 

0.89044 

1.14935 

(0.0000) 

(0.1572) 

(0.2176) 

(0.2884) 

(0.3575) 

(0.4387) 

(0.5233) 

(0.6174) 

(0.7478)* 

Depth- 
expanded, 
eq. (48) 

0.00000 

0.01682 

0.06861 

0.15723 

0.28417 

0.45001 

0.65358 

0.89041 

1.14929 

(0.0000) 

(0.1572) 

(0.2176) 

(0.2884) 

(0.3575) 

(0.4395) 

(0.5242) 

(0.6200) 

(0.7530) 


* Most accurate prediction (closest to SPAR). 


Table 2(d). Comparisons of deflections (in inches) calculated from SPAR and from different displacement 
transfer functions; tapered cantilever tubular beam; c n /c Q = 1/4; n = 16 . 


Theories 

yo 

Y2 

Y4 

Y6 

ys 

yio 

yi2 

yi4 

yi6 

(% error) 

SPAR 

(reference) 

0.00000 

(0.0000) 

0.01912 

(0.0000) 

0.07598 

(0.0000) 

0.17923 

(0.0000) 

0.33816 

(0.0000) 

0.56506 

(0.0000) 

0.87454 

(0.0000) 

1.28196 

(0.0000) 

1.78417 

(0.0000) 


Piecewise linear strain representations 


Nonuniform, 
eq. (4) 

0.00000 

(0.0000) 

0.01739 

(0.0970) 

0.07364 

(0.1312) 

0.17612 

(0.1743) 

0.33427 

(0.2180) 

0.56011 

(0.2774) 

0.86835 

(0.3469) 

1.27375 

(0.4602) 

1.77144 

(0.7135) 

Second- 
order, 
eq. (6) 

0.00000 

(0.0000) 

0.01739 

(0.0970) 

0.07364 

(0.1312) 

0.17611 

(0.1749) 

0.33426 

(0.2186) 

0.56009 

(0.2786) 

0.86829 

(0.3503) 

1.27365 

(0.4658) 

1.77124 

(0.7247) 


Piecewise nonlinear strain representations 


Improved, 
eq. (27)* 

0.00000 

(0.0000) 

0.01738 

(0.0976) 

0.07362 

(0.1323) 

0.17609 

(0.1760) 

0.33424 

(0.2197) 

0.56010 

(0.2780) 

0.86850 

(0.3385) 

1.27462 

(0.4114) 

1.77473 

(0.5291)* 

Log- 

expanded, 
eq. (33) 

0.00000 

(0.0000) 

0.01738 

(0.0975) 

0.07362 

(0.1323) 

0.17608 

(0.1766) 

0.33422 

(0.2208) 

0.56007 

(0.2797) 

0.86845 

(0.3413) 

1.27452 

(0.4170) 

1.77468 

(0.5319) 

Depth- 
expanded, 
eq. (48) 

0.00000 

(0.0000) 

0.01738 

(0.0975) 

0.07362 

(0.1323) 

0.17608 

(0.1766) 

0.33422 

(0.2208) 

0.56008 

(0.2791) 

0.86845 

(0.3413) 

1.27453 

(0.4164) 

1.77458 

(0.5375) 


* Most accurate prediction (closest to SPAR). 
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Table 2(e). Comparisons of deflections (in inches) calculated from SPAR and from different 
Displacement Transfer Functions; tapered cantilever tubular beam; c n /c 0 = 0.5/4; n = 16 . 


Theories 

Yo 

Y2 

Y4 

Y6 

ys 

YlO 

y 12 

yi4 

yi6 

(% error) 

SPAR 

0.00000 

0.01938 

0.07872 

0.19047 

0.37082 

0.64535 

1.05517 

1.67292 

2.59712 

(reference) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

Piecewise linear strain representations 

Nonuniform, 

0.00000 

0.01770 

0.07649 

0.18752 

0.36721 

0.64077 

1.04963 

1.66517 

2.57675 

eq. (4) 

(0.0000) 

(0.0620) 

(0.0859) 

(0.1136) 

(0.1390) 

(0.1763) 

(0.2133) 

(0.2984) 

(0.7843) 

Second- 

order, 

eq. (6) 

0.00000 

0.01770 

0.07648 

0.18751 

0.36719 

0.64071 

1.04949 

1.66478 

2.57544 

(0.0000) 

(0.0620) 

(0.0862) 

(0.1140) 

(0.1398) 

(0.1787) 

(0.2187) 

(0.3134) 

(0.8348) 

Piecewise nonlinear strain representations 

Improved, 

0.00000 

0.01768 

0.07644 

0.18741 

0.36700 

0.64038 

1.04905 

1.66527 

2.58806 

eq. (27)* 

(0.0000) 

(0.0655) 

(0.0878) 

(0.1178) 

(0.1417) 

(0.1914) 

(0.2356) 

(0.2946) 

( 0 . 3488 )* 

Log- 

expanded, 
eq. (33) 

0.00000 

0.01769 

0.07644 

0.18740 

0.36697 

0.64030 

1.04887 

1.66475 

2.58595 

(0.0000) 

(0.0651) 

(0.0878) 

(0.1182) 

(0.1482) 

(0.1944) 

(0.2426) 

(0.3146) 

(0.4301) 

Depth- 
expanded, 
eq. (48) 

0.00000 

0.01768 

0.07644 

0.18740 

0.36698 

0.64032 

1.04891 

1.66489 

2.58719 

(0.0000) 

(0.0655) 

(0.0878) 

(0.1182) 

(0.1479) 

(0.1937) 

(0.2410) 

(0.3092) 

(0.3823) 


* Most accurate prediction (closest to SPAR). 


Table 2(f). Comparisons of deflections (in inches) calculated from SPAR and from different 
Displacement Transfer Functions; tapered cantilever tubular beam; c n /c 0 = 0.25/4; n = 16 . 


Theories 

yo 

Y2 

Y4 

Y6 

ys 

yio 

yi2 

yi4 

yi6 

(% error) 

SPAR 

(Reference) 

0.00000 

(0.0000) 

0.01952 

(0.0000) 

0.08017 

(0.0000) 

0.19668 

(0.0000) 

0.38984 

(0.0000) 

0.69566 

(0.0000) 

1.18157 

(0.0000) 

2.00286 

(0.0000) 

3.58289 

(0.0000) 


Piecewise linear strain representations 


Nonuniform, 
Eq. (4) 

0.00000 

(0.0000) 

0.01786 

(0.0463) 

0.07802 

(0.0600) 

0.19387 

(0.0784) 

0.38651 

(0.0929) 

0.69143 

(0.1181) 

1.17563 

(0.1658) 

1.99425 

(0.2403) 

3.53790 

(1.2557) 

Second- 
order, 
eq. (6) 

0.00000 

(0.0000) 

0.01786 

(0.00463) 

0.07801 

(0.0603) 

0.19385 

(0.0790) 

0.38648 

(0.0938) 

0.69135 

(0.1203) 

1.17539 

(0.1725) 

1.99331 

(0.2665) 

3.53116 

(1.4398) 


Piecewise nonlinear strain representations 


Improved, 
eq. (27)* 

0.00000 

(0.0000) 

0.01784 

(0.0469) 

0.07795 

(0.0620) 

0.19369 

(0.0835) 

0.38615 

(0.1030) 

0.69068 

(0.1390) 

1.17400 

(0.2113) 

1.99125 

(0.3240) 

3.57152 

( 0 . 3173 )* 

Log- 

expanded, 
eq. (33) 

0.00000 

(0.0000) 

0.01784 

(0.0469) 

0.07794 

(0.0622) 

0.19368 

(0.0837) 

0.38609 

(0.1047) 

0.69057 

(0.1421) 

1.17372 

(0.2191) 

1.99044 

(0.3466) 

3.56071 

(0.6191) 

Depth- 
expanded, 
eq. (48) 

0.00000 

(0.0000) 

0.01784 

(0.0469) 

0.07794 

(0.0622) 

0.19368 

(0.0837) 

0.38611 

(0.1041) 

0.69058 

(0.1418) 

1.17373 

(0.2188) 

1.99026 

(0.3517) 

3.56740 

(0.4323) 


* Most accurate prediction (closest to SPAR). 


The deflections calculated from equations (4), (6), (27), (33), and (48), and listed in tables 2(a)-2(f) 
are plotted respectively in figures 9-13 (solid curves with solid circular symbols) for visual comparisons 
with the corresponding SPAR deflection curves (dashed curves with open circular symbols). Note that by 
using the piecewise nonlinear strain representations (figs. 11-13) instead of the piecewise linear strain 
representations (figs. 9 and 10), the shape-prediction accuracies, especially for the smaller depth ratios, 
c n /c 0 = {0.25/4, 0.5/4}, were greatly improved. Pictorially, the deflection curves calculated from 
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equations (27), (33), and (48), and shown respectively in figures 11-13, fell extremely close to the 
corresponding SPAR deflection curves because the piecewise nonlinear strain representations were used 
[see the actual values of the prediction errors indicated within the parentheses in tables 2(a)-2(f)]. 

Shape-Prediction Errors 

In the shape-prediction error studies, the SPAR-generated deflections were used as the reference 
yardsticks with which to measure the shape-prediction errors. The shape-prediction error studies focused 
only on the beam tips because of peak deflections. The beam-tip shape-prediction error, 
A y n (=| y n -yf A *\) (in.) for each type of beam is defined as the magnitude of difference between the 

SPAR-generated beam-tip deflection, y f *** , and the corresponding beam-tip deflection, y n , calculated 

from equations (4), (6), (27), (33), and (48). From data listed in Table 2(a)-2(f), the calculated beam-tip 
shape-prediction errors, A y n , are compared in table 3. 


Table 3. Comparisons of beam-tip shape-prediction errors, A_y n (in.) of different Displacement Transfer 
Functions; tapered cantilever tubular beams with depth ratios 4/4 >c n /c Q >0.25/4. 


C J C 0 

V®" 

in. 


Beam-tip shape-prediction 

error, A y n , in. 


Piecewise linear 


Piecewise nonlinear 

Nonuniform, 
eq. (4) 

Second - 
order, 

eq. (6) 

Improved, 
eq. (27) 

Log- 

expanded, 
eq. (33) 

Depth-expanded, 
eq. (48) 

4/4 

0.71162 

— 

0.00919 

— 

0.00913 

0.00913 

3/4 

0.87634 

0.00923 

0.00922 

0.00880 

0.00879* 

0.00880 

2/4 

1.15801 

0.00973 

0.00975 

0.00871 

0.00866* 

0.00872 

1/4 

1.78417 

0.01273 

0.01293 

0.00944* 

0.00949 

0.00959 

0.5/4 

2.59712 

0.02037 

0.02168 

0.00906* 

0.01117 

0.00993 

0.25/4 

3.68289 

0.04499 

0.05173 

0.01137* 

0.02218 

0.01549 


* Most accurate. 


The beam-tip shape-prediction errors, A y n , listed in table 3 are plotted in figure 14 for visual 

comparisons of the shape-prediction accuracies of the Displacement Transfer Functions (4), (6), (27), 
(33), and (48). One observes from figure 14 that, for the uniform beam (c n l c 0 = 4/4), changing from 

piecewise linear to piecewise nonlinear strain representations has no effect in reducing the shape- 
prediction errors because it is immaterial to use either piecewise linear or piecewise nonlinear function to 
represent the linear (straight-line) strain curve for the uniform beam (figs. 5 and 8). Note also that the 
values of the beam-tip shape-prediction errors, Ay n , are practically insensitive to the change of the depth 

ratio, c n /c 0 , and stayed almost unchanged in the region, 4/4 >c n /c 0 > 2/4, but started to increase 
rapidly as the depth ratio, c n /c 0 , decreases, and reached the maximum values of 0.04499, 0.05173, 
0.01137, 0.02218, and 0.01549 in. at the minimum depth ratio, c n /c 0 = 0.25/4, respectively for equations 
(4), (6), (27), (33), and (48) (table 3). Because at lower depth ratios, c n lc 0 < 1/4, the strain curves are 

highly nonlinear (fig. 8), and using the piecewise nonlinear strain representations [equations (27), (33), 
and (48)], instead of using the piecewise linear strain representations [equations (4) and (6)], greatly 
reduced the beam-tip shape-prediction errors (fig. 14). Note also in figure 14 that all the error curves bent 
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upward (increasing errors) when the depth ratio decreased toward the minimum depth ratio, 
c n /c 0 = 0.25/4. This upward bend is most severe for the piecewise linear strain representations 
[equations (4) and (6)], and to a lesser degree for the piecewise nonlinear strain representations [equations 
(27), (33), and (48)]. The increasing errors at low depth ratios, c n /c 0 = {0.25/4, 0.5/4}, can be attributed 

to the relatively poor strain representations of the sharply-bent segments of the strain curves (figs. 5 and 
8). The beam-tip shape-prediction errors plotted in figure 14, however, are very small and under 0.06 
inches. The Improved Displacement Transfer Function (27) has the smallest error (0.01 137 in.). 

Normalized Shape-Prediction Errors 

The beam-tip shape-prediction errors, A y n , of different Displacement Transfer Functions can be 
better compared if each A y n is normalized by the corresponding SPAR-generated beam-tip deflection, 
yf >AR , (prediction-error measuring yardsticks). The normalized beam-tip shape-prediction errors, 
A y n , of different Displacement Transfer Functions are compared in table 4. 

Table 4. Comparisons of normalized beam-tip shape-prediction errors (A^/y^^xlOO 0 /)) of different 
Displacement Transfer Functions; tapered cantilever tubular beams with depth ratios, 4/4 > c n /c Q > 1/4. 



Normalized beam-tip shape-prediction 

error, A yJy s „ PAt , 

XI 00% 


Piecewise linear 


Piecewise nonlinear 

C J C 0 

Nonuniform, 

Second- 

order, 

Improved, 

Log- 

expanded, 

Depth-expanded, 


eq. (4) 

eq. (6) 

eq. (27) 

eq. (33) 

eq. (48) 

4/4 

— 

1.2914 

— 

1.2830 

1.2830 

3/4 

1.0532 

1.0521 

1.0042 

1.0030* 

1.0042 

2/4 

0.8402 

0.8420 

0.7522 

0.7478* 

0.7530 

1/4 

0.7135 

0.7147 

0.5291* 

0.5319 

0.5375 

0.5/4 

0.7843 

0.8248 

0.3488* 

0.4301 

0.3823 

0.25/4 

1.2557 

1.4398 

0.3173* 

0.6191 

0.4323 


*Most accurate. 


The normalized beam-tip shape-prediction errors, A y n / y^ PAl (in percent), listed in table 4 are plotted 
in figure 15 as functions of depth ratio, c n i c 0 , for visual comparisons of the shape prediction accuracies 

of different Displacement Transfer Functions. Note from figure 15 that the values of Ay n !y^ PAl 

associated with all types of Displacement Transfer Functions decreased almost linearly with decreasing 
depth ratio, c n /c 0 , and then [except for the Improved Displacement Transfer Function (27)] increased 

rapidly with further decreasing depth ratio, c n /c 0 (see fig. 14). The reason for the values of Ay n /y^ PAl 
decreasing with decreasing cj c Q in the region, 4/4 >c n /c 0 >1/4, is dividing the almost-unchanging 

values of A y n (see fig. 14) by the associated values of , which increase with decreasing depth ratio, 
cj c o (see column 2 of table 3). 
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At the low depth ratios, c n /c 0 ={0.5/4, 0.25/4}, the Improved Displacement Transfer Function (27) 
gave the lowest value of A y n / y^ PAl among all of the five Displacement Transfer Functions that were 
considered (fig. 15). At the lowest depth ratio, c n /c 0 = 0.25/4, the beam-tip shape-prediction errors of the 

Improved, Depth-expanded, and Log-expanded Displacement Transfer Functions (27), (48), and (33) are 
respectively 25%, 34%, and 49% of the error produced by the Nonuniform Displacement Transfer 
Function (4) (the most accurate equation, based on piecewise linear strain representations, see table 4). 
This is equivalent to error reductions of 75%, 66%, and 51%, respectively, by using the Displacement 
Transfer Functions (27), (48), and (33). Note from figure 15 that in the depth ratio region, 
4/4 >c n /c 0 > 1/4, all three of the Displacement Transfer Functions (27), (48), and (33) exhibited nearly 

equal shape-prediction accuracies. Similarly, in the same depth ratio region, 4/4 >c n /c Q > 1/4, the error 

curves of the Displacement Transfer Functions (4) and (6) (based on piecewise linear strain 
representations), also collapsed into a single curve. Overall, the prediction errors of the five Displacement 
Transfer Functions are very small: only in the ranges of (0.7135-1.4398)% based on piecewise linear 
strain representations and (0.3173-1.2830)% based on piecewise nonlinear strain representations. 

DISCUSSION 

For shape predictions of complex aircraft wings (which are constructed with stiffened panels) using 
the Displacement Transfer Functions, one must first locate the local cross-sectional neutral axis using 
pairs of lower and upper surface strains, and then determine the local depth factors, c f (i = 1 , 2 , 3 ,...., ri) 

(see Appendix I), which are needed for the deflection calculations. The surface strain values will change 
with the material or structural properties. The above approach was applied in the shape predictions of an 
Ikhana wing (ref. 3) and a Global Observer wing (ref. 18), for which the local neutral axes are not at the 
half-depth of the wing cross sections. 

Concurrently with the writing of this technical paper, large-scale ground load tests of the Global 
Observer (175 -ft wingspan) were carried out at the NASA Dry den Flight Loads Laboratory (ref. 18). The 
measured strain test data were later used in the previously listed Displacement Transfer Function (5), 
(suitable for the slightly-tapered Global Observer wing) to calculate the wing deflections. The calculated 
values were then compared with the real-time photogrammetry data as shown in figures 16 and 17 (taken 
from ref. 18). Note that equation (5) proved to be very accurate, having prediction errors in the range of 
(0.09-1.72)% at the wingtip, which was at nearly 14 ft of deflection. 

CONCLUDING REMARKS 

New Improved Displacement Transfer Functions were formulated using the piecewise nonlinear 
representations of surface strain distributions. Four different Displacement Transfer Functions were 
formulated for both uniform and nonuniform beams: 

• Improved Displacement Transfer Function for uniform beams 

• Improved Displacement Transfer Function for nonuniform beams 

• Log-expanded Displacement Transfer Function for nonuniform beams 

• Depth-expanded Displacement Transfer Function for nonuniform beams. 

The key findings resulting from using the new Displacement Transfer Functions to improve the 
accuracies of structure shape predictions of tapered cantilever tubular beams are: 

1. By using the piecewise nonlinear strain representations, beam-tip prediction errors could be 
greatly reduced by up to 75% for the peak taper angle case. 
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2. For smaller taper angles, changing the strain representations from piecewise linear to piecewise 
nonlinear improved the shape-prediction accuracies along the entire beam. 

3. For larger taper angles, changing the strain representations from piecewise linear to piecewise 
nonlinear improved the prediction accuracies only in the beam tip regions where the strain 
gradients change signs, but did not improve the prediction accuracies in the inboard regions 
where the strain gradients change very slowly. 

4. For larger taper angles, among the three new Displacement Transfer Functions formulated for 
nonuniform cantilever beams, the Improved Displacement Transfer Function provided the most 
accurate shape predictions. 

5. For smaller taper angles, the shape-prediction accuracies of the Improved, Depth-expanded, and 
Log-expanded Displacement Transfer Functions are practically the same. 

6. For the uniform beams, there is no benefit in the prediction error reductions by changing from 
linear strain representations to nonlinear strain representations. As the depth ratio decreases (that 
is, increasing taper angles), however, the benefits of prediction error reductions gradually 
increased and finally reached the peaks at the minimum depth ratio (the highest taper angle). 
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FIGURES 



Figure 1. Small segment of deformed uniform embedded beam for relocating local curvature R(x) to 
associated surface bending strain s(x). 



Figure 2. Embedded beam (cross section of any structure, including tubular beam) along the 
strain-sensing line with evenly distributed strain-sensing stations. 
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1 = 100.5 in. 
t =0.02296 in. 
c Q = 4 in. 

c n = 4, 3, 2, 1, 0.5, 0.25 in. 


Nodes 3673 

Four-node elements 3636 
Three-node elements 36 


Figure 3. Finite-element model for tapered cantilever tubular beam subjected to tip vertical load 
P= 1001b. 
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Strain-sensing stations, i 6.5 

3 ’ 120159 

Figure 4. SPAR-generated surface strain data points fitted with piecewise linear strain curves calculated 
from equation (9) for different tapered cantilever tubular beams discretized with n = 8 domains; 
P= 1001b. 



Figure 5. SPAR-generated surface strain data points fitted with piecewise linear strain curves calculated 
from equation (9) for different tapered cantilever tubular beams discretized with n = 16 domains; 


P= 100 lb. 
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c ( t \_ c _ 3E /-1 ~ 4E f + E /+1 ( r - y . \ . e /-1 “ 2£ /‘ + e /+1 fr _ T . \2 

( } ~ '- 1 2A7 ( m) + 2(A/)2 — ( m) 


Figure 6. Local second-order curve fitting of three adjacent strain points. 



120162 


Figure 7. Extrapolation of three adjacent strains {£ n _ 2 , £ n- i> £ n} t° obtain fictitious strain, £ n+1 . 
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Strain-sensing stations, i 14.5 

11201163 

Figure 8. SPAR-generated surface strain data points fitted with piecewise linear strain curves calculated 
from equation (10) for different tapered cantilever tubular beams discretized with n = 16 domains; 
P= 1001b. 
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Strain-sensing stations, / 

120164 

Figure 9. Comparisons of deflections calculated from SPAR program and from Nonuniform 
Displacement Transfer Function [equation (4)] (based on piecewise linear strain representations) for 
tapered cantilever tubular beams; P = 100 lb. 



Strain-sensing stations, / 

120165 

Figure 10. Comparisons of deflections calculated from SPAR program and from Second-order 
Displacement Transfer Function [equation (6)] (based on piecewise linear strain representations) for 
tapered cantilever tubular beams; P= 100 lb. 
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Figure 11. Comparisons of deflections calculated from SPAR program and from Improved Displacement 
Transfer Function [equation (27)] (based on piecewise nonlinear strain representations) for tapered 
cantilever tubular beams; P = 100 lb. 



Strain-sensing stations, / 

120167 

Figure 12. Comparisons of deflections calculated from SPAR program and from Log-expanded 
Displacement Transfer Function [equation (33)] (based on piecewise nonlinear strain representations) for 
tapered cantilever tubular beams; P= 100 lb. 
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Figure 13. Comparisons of deflections calculated from SPAR program and from Depth-expanded 
Displacement Transfer Function [equation (48)] (based on piecewise nonlinear strain representations) for 
tapered cantilever tubular beams; P= 100 lb. 



Depth ratio, c n /c 0 

120169 

Figure 14. Comparisons of beam-tip deflection prediction errors of different Displacement Transfer 
Functions based on piecewise linear and piecewise nonlinear strain representations. 
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Depth ratio, c n /c 0 

120170 

Figure 15. Comparisons of percent beam-tip deflection prediction errors of different Displacement 
Transfer Functions based on piecewise linear and piecewise nonlinear strain representations. 
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Figure 16. Comparison of predicted and measured wing deflections along No. 3 strain-sensing line for 
different loading levels (ref. 18). 



Figure 17. Comparison of predicted and measured wing deflections along No. 4 strain-sensing line for 
different loading levels (ref. 18). 
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APPENDIX A 


DERIVATIONS OF IMPROVED SLOPE AND DEFLECTION 
EQUATIONS FOR UNIFORM BEAMS 

Appendix A presents the details of integrations of the slope equation (12) and the deflection equation 
(16) for the uniform beams to obtain the final mathematical forms given respectively by equations (14) 
and (18). 

Slope Equations 

The slope, tan 0(x ) □ of the uniform beam in the small domain, x ;1 <x<x f , between the two 
adjacent strain-sensing stations, {jy , , x, } , can be obtained by carrying out the integration of equation 
(12), which is reproduced below: 

taa0(x) = x ^i^=\ X ^^- dx + tan 0 t _ x ; (x, , <x<x,) (A . ]} 

v -v- 7 Slope at x f _j ' — ^ — * 

Slope increment Eq. (2) 


The nonlinear distribution of strain, s{x ), in the domain, x i X <jc<jc-, between the two adjacent 
strain-sensing stations, { x i _ l ,x i }, described by equation (10), is duplicated below: 


£ (x) = - 


3£T , - 4£T +s m , . £ t , -2t- + e . 2 

11 1 i+1 (x-x_ 1 )+ 11 i+1 (x-x_,) 2 (A-2) 


2M 


2(A/) 


(*/_i < x < x,) 


In view of equation (A-2), equation (A-l) can be integrated to yield (ref. 17): 


tan 0(x) = 1 f 

C Jx n 


£ i-! ~ 


3 £ t i 4g t - + g /+ ] \ . £ l i 2i- L + £ [iX ( _ -\ 2 

2 A/ 1 ilJ 2(A If 


L (x-x i _,) 2 


ftv + tanft 


i-l 


£^(x-x h )- 


4A/ 


-(x-x i _ i f+- 


6(A If 


+ tan ft 


i-l 


3 £ i - 1 4f i + £ m ^ , £ i-\ 2 -£i + £ M( x _ x _ ys 

(x,_, <X<X,) (A-3) 


which is equation (13) in the text. At the strain-sensing station x f , one can write x — x, , = A/, and 
equation (A-3) yields the slope, tan ft [ = tanft(x ( )] , at the strain-sensing station. x i as: 


„ A l 

tan ft = — 


3 £ !-i - 4i", + £ l+l , £ i _ l -2£ l + £ i+l 

A-l . ' r 


+ tan ft t ; (/ = 1, 2, 3,. ...,«) (A-4) 


After grouping terms, equation (A-4) takes on the following form: 
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(A-5) 


which is equation (14) in the text, the Improved slope equation for the uniform beams. 

Deflection Equations 

The deflection, y(x ), of the uniform-depth beam in the small domain, x zl <x<x r between the two 
adjacent strain-sensing stations, { x il ,x i }, can be obtained by integrating the slope equation (A-3) 

[or equation (13)] with constant of integration determined by enforcing the continuity of deflection at the 
inboard adjacent strain-sensing station, jc- 1? as: 


in which y iX is the deflection at the strain-sensing station x i X . 

Substituting equation (A-3) into equation (A-6), and carrying out integrations, one obtains (ref 17): 


which is equation (17) in the text. At the strain-sensing station, x v one can write x t ~x i X = A/, and 
equation (A-7) gives the deflection, yf= at the strain-sensing station, x 2 , as: 


X-^) = fAan^)^+ ; (x ;1 <x<x,) 


(A-6) 


Eq. (A-3) Deflection 


Deflection 






+X-i+( A/ ) tan ^-i (A' 8 ) 

(* = 1,2,3,.— ,w) 


After grouping terms, equation (A-8) takes on the form: 



(A-9) 


which is equation (18) in the text, the Improved deflection equation for the uniform beams. 
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APPENDIX B 


DERIVATIONS OF IMPROVED DISPLACEMENT TRANSFER 
FUNCTION FOR UNIFORM BEAMS 

Appendix B presents the derivation of the final dual summation form of the deflection equation for 
the uniform cantilever beams. Equations (14) [or equation (A-5)] and (18) [or equation (A-9)] are 
duplicated below as equations (B-l) and (B-2), respectively, and can be combined into one equation with 
two summation terms as derived below. 

tan 0 i = A- (5/;, , + 8c, - c w ) + tan () i , ; (/ = 1,2,3,...., n) (B-l) 

12c 

X=^^( 7f i a + 6£ ;-%i)+X-i+ A/tan ^-i ; (/= 1,2, 3,. (B-2) 

In view of equation (B-l), equation (B-2) may be written for different indices as: 


For i = 1 : 


» = (7 s o + 6c, - c 2 ) + y 0 + Al tan<9 0 

24c 


(B-3) 


For / = 2 : 


y 2 = ^-(7c, + 6c 2 - c 3 )+y + A/tan 6> 

24 c 

(A If (Al) 2 

= — ~ ~ (7* + 6c 2 - c 3 ) + — — (7c 0 + 6c, - c 2 ) + y 0 + Al tan 0 O 
24 c 24c 


+ ^T“ ( 5 ^o + 8c, - c 2 ) + A/ tan <9 0 
12c 

V V ' 

Al tan Gy 


(B-4) 


After grouping terms, equation (B-4) becomes: 

y, = (7£, +6^ (7 £ , +6£, -£ 2 )+(2-l)l^f (5s 0 +8£, -£,) 

24c 24c 12c (B-5) 

+y 0 +2A/tan6> 0 
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For i = 3: 


y 3 =^-(l£ 2 +6s i -s A )+y 2 + A/tan 0 2 
24 c 


(A If 

24c 


( 1 s 2 + 6 s 3 - s a ) 


(A If (A/) 2 (A/) 2 , „ 

+— — (76r, +6s 2 -£ 3 )+— — (7s 0 +6£r, — £ 2 )+(2—l)— — (5£r 0 +8cj -£ 2 )+y 0 +2A/tan 6> 0 
24c 24c 12c 


^2 


+-— — — (5^ +8g 2 -g 3 ) + ^^ (5g 0 +8g, — g 2 ) + A/tan6> 0 
12c 12c 


(B-6) 


A/ tan ^ 


After grouping terms, equation (B-6) becomes: 


(A/) 2 (A/) 2 (A/) 2 

y 3 = 2^ (2^2 ^^3 — ^4) ( 7f i ^^2 — ^3) 4 " + 66 1 ] — £ 2 ) 

+ (3 - 2) ^-(5c, +S£ 2 -£ 3 )+(3- 1) (5f 0 + 8c, - £ 2 ) + y 0 + 3 A/ tan 0 O 


(B-7) 


Observing the indicial behavior, equation (B-7) can be generalized for index /, and the deflection, j/_, 
can be expressed in dual summation form as: 


(\ A 2 1 (\f \ 2 11 

y< = +6e j - £ j + i) + yo +^rL<f-jX 5€ j-i +8 s -%,)+(OVtan6> 0 (B8) 


Contributions from deflection terms 


Contributions from slope terms 

(* -1,2,3,-— ,«) 


which is equation (19) in the text, the Improved Displacement Transfer Function for the uniform beams. 
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APPENDIX C 


DERIVATIONS OF IMPROVED SLOPE AND DEFLECTION 
EQUATIONS FOR NONUNIFORM BEAMS 

Appendix C presents the details of integrations of the slope equation (20) and the deflection equation 
(24) for the nonuniform beams to obtain the final mathematical forms given respectively by equations 
(22) and (26). 


Improved Slope Equations 

The slope, tan£>(x)D of the nonuniform beam in the small domain, x ; , <x<x r between the two 
adjacent strain-sensing stations, j x ( , , x ; | . can be obtained by carrying out the integration of equation 
(20) which is reproduced below as equation (C-l): 


tan 0{x) = [ ^-^-dx + tan Q t _ x = f ^^t&+tan^_, ; (x..<x<x) 

Jx. . nv % / Jx. , ft rl 


£(x) 


,x <-i cbc 

' Tv ' 

Slope increment 


Slope at Xj 


,x t - 1 c{x) 

Eq (3) 


(C-l) 


The distributions of beam depth factor, c(x) [equation (8)], and bending strain, x(x) [equation (10)], in 
the small domain, x ; , <x<x r between the two adjacent strain-sensing stations, \ x t , , x | , are 
respectively duplicated below as equations (C-2) and (C-3): 


c( x ) = c l -i + (c i -c l _ 1 ) 11 ; (x, ,<*<*,) 


(C-2) 

s(x) = £ 3e ‘^~ 4£ ! +£ M( x X )+ £ ‘-l ~ 2£ i +£ M, x x f . 

W 11 2A1 { , l) 2(A If 1 , l) ’ 

(x, i <X<X,) 

(C-3) 


Substituting equations (C-2) and (C-3) into equation (C-l) yields: 


tan 0(x) = f 

1 


'is, . - 4s + s M . s . , - 2 s i + s M , . 2 

^ i A77 ( x ~ x > ■ ) + ‘ 2 


2A/ 


2(A/> 


c. , + 


2AI 


< x ~ x i- 1) 


dx + tenO i , (C-4) 


(A, — x < Xj ) 


Let 
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4 _ 3 ^,- 4 ^ + ^, 

(C-5) 

2A7 

^ _ £ i-i ~ 4- f i+1 

2(A/) 2 

(C-6) 

^ c — c , 


c = * '- 1 

(C-7) 

A/ 

£ = (x-x, j) 

(C-8) 


then equation (C-4) takes on the following simplified form: 


tan#(x) = dZ + tmO^ (C-9) 


After carrying out integration of equation (C-9), one obtains (ref. 17): 


tan = 


^ 1 °g e ( c /-i + c '#) + aI-A^c» + C€) 

v v 7 v v 7 

Integration of 1st term in eq. (C-4) 


~l4 


B 


Integration of 2nd term in eq. (C-4) 

\2 ~ , 2 


+ Cff - 2c,_,(c,_, + CZ) + c;_, log.(c,_, + CD 


Integration of 3rd term in eq (C-4) 

1 B 


+ tan# 


B 


+ B ' 1 I log K (c , +C£) + -Z+ (c^ + C£) - 2 —c l ,(c , + C0 


C 2 C 


2 C 3 


C 



. 1 B 2 

log,^ , 


„ 5 2 

+2 e c ‘ 


au^., 


+ 


+ 


H- 


A[c 2 c h - ACc i _ i + Be’ , ][log,(c,_, + CD - log, c„ ] + 
j + <"D' -<£,]- 2 Ac, , [(c„ + CD - c,_, ] + tan f ) : , 

A[ C ’c m - /ICc,., +«4,][log,(c,. 1 + CD-log,c 1 _,]+ ^ 
^ £ ( 2 <^ + C#” ) - 2 + tan 0 ,_, 

3j(CXi ~ ACc i-t +/fei,][l°g,(c,_| +CD-log,c,_,]+^ 
iA[_2c,^ + CD] + tane M 


(C-10) 
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In view of the definitions from equations (C-5), (C-6), (C-7), and (C-8), equation (C-10) takes on the 
following form: 


tan0(x) = 


(A/)’ 




( C i C i- 1) , (? £ i-X ^ £ i ~*~ £ i+l) ( C i C i 1 ) „ , (^7-1 2f r + f i+ i) 


logj C M + 


(A/) 2 2 A/ 

-(x-x, | )]-]°g f ,C i _ 1 


A/ 


A/ C,t + 2(A/) 2 

(3^-_i-4£r+£r i+1 ) A/ 


2A/ 


(x-x, J 


+ 


1 (^,-2^+^) (A If 


2 2(A/) 2 (c,-^) 2 


-2c, j(x - x M ) + - x t _ x ) 2 


A/ 


+tan# 


A/ 


, v , £,_i (cr c /-i) 2 + h 3 ^ i - 4s i + %i)( c , -C mX-i + h*', , - 2s, + s l+l )cf , 
(c -c, ,) L 2 2 


lo g e | fi-i + Cf i C ‘ ' (* ~ JCf-i ) I ~ log. c,_ t 


— 4 -I- /* 

/-t / V(x-x i i )+" 11 ‘ 


Al 

£ , — 2c + £ 


2{c-c i _ l ) 


4 ( c i L 


-2c i . 1 (x-x i _ l)+ ^(x-x i _ 1 ) 2 

Al 


+ tant? 


(x M <x<x,) (C-ll) 


In equation (C-l 1), the two terms, { I l ,J 2 }, can be simplified through grouping terms as follows: 


2 1 1 2 

/j = £ 1} (C' ~c i } ) + —Q £ i~i ~ 4 £ , + e M ){c t -c t _ x ) Ci _, +-( £ i-i -2 e, + £ M )c t i 


( C i C 1 l) 1 2^ C ‘ C ' ^ C i -1 + 2 C ' _1 


1 [2(c z C, l) 6 ', 1 + C f_l + [( C , C , l) 6 ’;-l + C /-1 J^I+1 


2 ~ 2 3 3 2 1 2 1 s 1 

= yc, - 2c j c i _ l + c M + -cp^ - -c l X + -c l X J c,, - (2c, - c,_, )c, x s, + -cp t _ x c M 

= ' [( 2 C - C ,-l) C , £ ,-l - 2 ( 2c , - Ci-iPh*, + C i C il £ M ] 


(C-l 2) 


/, = _ (^ -^i+0 (x _^ i)+ ^-2s, + s M 


2(c, -c H ) 

1 


4(c i -c,_ 1 ) 2 


(c — c t ) 9 

-2c H (x-aO+ ^(x-x,.,) 2 
Al 


2 (c,-c M ) 2 


[(- 3 ^-i + 4 C - %i )(C - c ,-, ) - (Vi - 2«i + K i ] (* - ) 


+ 


c , -2c +c + , (x-x, t ) 2 
4(6- -C,,) Al 
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= zr — 1 — ;^{[ _3 (e -C-i)-C-iKi + [ 4 (e - C-.) + 2c ,. ] C + Re -c f - 1 )-c i _ 1 ]^ + i}(x-x I _i) 

2(C“C-i) 

_ L ff i _,-2«i + «^ (x-x, ,) 2 


4(C-C_i) 

1 

2(C _ C-i) 2 


A/ 


[( 2 C-i -3c, )^, , + 2(2c i -c i _ l )e i -c i £ M ](x-x i _ l )+ 


C-i ~ 2gj + C + i (x-x,,,) 2 

4(c ; - C|1 ) A/ 


(C-13) 


Substitutions of equations (C-12) and (C-13) into equation (C-l 1) yields: 


M 


hmO(x)= [(2c -2 c_ 1 ^)+c£_ 1 C + i] 

2(C-C- 1 ) 


lo 8 e I C l + C ‘ > (X Xj_!> I - log, c,_, 


A/ 


— 1 2 [(2c,_ 1 -3c K_, + 2(2c 

2(C-C-i) 


+ *m~ 2«, +«„(*-*„) +|ang 


4(c,-c,.,) 


A/ 


(Xj_, <x<x,) (C-l 4) 


At the strain-sensing station, x r one can write (x j - x^ t ) = Al , and equation (C-l 4) gives the slope, 
tan 0\ = tan 0(x i ) | , at the strain-sensing station, x r as 


A/ r i c 

ton °f = yT (2C - XCC-i - 2C ,C>+ CC-iCu J log, — 

2(c-c,-_ 1 ) C-i 


A, — , 2 [( 2 C i -3c)C i +2(2c, -C_i)C -CC + i] + £ ' * — + ^ +1 + tan 6> , 
z vC Ci) 'HC Ci) 


(C-l 5) 


0 = 1, 2, 3,. 


The term I 3 in equation (C-15) can be simplified as follows: 




A/ 

2(C“C iX 


[( 2 c_i - 3 e)e i + 2 ( 2 c, - e i)e - cc„ ] + 


2s ,+£,. 


4(c — C-i) 


A / 
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M 

2 (*,-*i-i) a 


Al 


2c, 1 -3c,+l(c,-c M ) 


1 


*1-1 + [ 2 ( 2C , “ C i X ) - (*l " C H )] *1 




[(5c, - 3c , )c , - 2(3c, - c_,)c + (c, + c,)c +1 ] 


(C-16) 


In view of equation (C-16), equation (C-15) takes on the final form as: 


A l c 

taT1 Oi = — TJ [(2 c t - C, , )(c^, , - 2c l l c,)+ c,c,_,c i+1 ] log, 

2(*,--*«) *« 

- — — , 2 [(5c f - 3*1-1 K, - 2(3*, - c «-i X + (*, + c «-i Ki ] + 1311 0 i-x 

4(c,-c,i) 

0 = 1,2, 3,. (C-17) 


which is equation (22) in the text, the Improved slope equation for the nonuniform beams. 


Improved Deflection Equations 

The deflection, y(x ), of the uniform beam in the small domain, x iX <x<x r between the two 
adjacent strain-sensing stations, { x il ,x i }, can be obtained by integrating the slope equation (C-14) with 

the constant of integration determined by enforcing the continuity of deflection at the inboard adjacent 
strain-sensing station, x t l , as 


y(x) = | tan0(x )dx + y iX ; (x^ <x< x t ) 

Eq. (C-14) Deflection 

atX i-4 


(C-18) 


Substituting equation (C-14) into equation (C-18), and carrying out the integrations (ref. 17), one obtains: 
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y(x ) . 


A/ 


2(c,-c w ) 


3 {[(2c, “ C i - lc i-\ £ i ) + c , c ,-i £ , + i] 


A/ 


A/ 


-(* - .T ; ) - c,_, logc w ) - (Jr - log, c„ 


(x - A,) 


4(c,-c h ) 


I [(2c M " 3 c,)£ m + 2(2 c, -c m )£, - c,£ (+1 ](a - a,.) 2 


g M ~2g, ■ + k~A,)- 

12(c t -c,.,) A/ 


+ (A-A,.)tan0 M + y._ 


(*, L < a < A ; ) (C-19) 


At the strain-sensing station, a ; , one can write x i —x [ , = A/ , and equation (C-19) gives the 
deflection, y,[= v ; (a, ) | . at the strain-sensing station, x r as follows: 

(A If fr , 

y, = — 7? ( 2 C - c ,-i - 2 c^s ,) + J 

2(C- C ,i> 

x [c log, c, - (c, - c H )- c f _! log c,_, - (c, - c,_! )log, c, , ] } 

(AO 2 


+- 


+ 


(A /) 2 

2 (c- Cl J 

(AO 2 


(2c,_i 3c ( )g,_i +2(2c ( C i )C g/gi+i + ^(g/-i 2^: + £ nl Xc C ,) 
+y i ,+A/tan6> ;1 

[(2c, - C-.XCC,-, - 2c m £,) + c M+l ][ Ci (log e c - log, c,)- (c - c M )] 




" ‘33 

+y ;1 + A/ tan , 

After grouping terms, equation (C-20) becomes: 


Cl - — 1 c, i + ( 4c, - 2c,_j - |c f + |ci_! ]c,+[-c, + --— k. 


, ■ 3 3 i-m 

0 = 1, 2, 3,. ...,«) (C-20) 


y = 


(AO 2 


2( c i 


-[(2c, -c i _ l )(c i e i _ l -2c i _ l £ i )+c i c i _ 1 e M ] 


c i lo 8 e — -( c , -c, i) 


(A /) 2 


[(8c, - 5c - 2(5c- - 2 c_ 1 )£ + (2c + c^s^] +y l _ l + M tantf ;1 


i2(<*-c M y 

0 = l»2,3,....,w) (C-21) 

which is equation (26) in the text, the Improved deflection equation for the nonuniform beams. 
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APPENDIX D 


DERIVATIONS OF IMPROVED DISPLACEMENT TRANSFER 
FUNCTION FOR NONUNIFORM BEAMS 


Appendix D presents the derivations of the Improved slope and deflection equations (22) [or equation 
(C-17)] and (26) [or equation (C-21)] [duplicated below as equations (D-l) and (D-2), respectively] that 
will be combined into one equation in dual summation form. 


tan6> = — -y[(2 c, -c H )(c^ H -2c^e l ) + cf M £ M \\og'-*- 


2( c , _c i-i) 

A/ 

'«c t -c„y 


u /-l 


[(5c - 3c, 2(3c, - c X + (c + e , )c , , ] + tan 6> , 

(/ = 1,2, 3,. 


(D-l) 


A,= 


(AO 2 


2(c i -c I _ 1 ) 4 

(AQ 2 


[(2c, c i_iX c A-i 2c, t c,) + c, c, , c, , , j 


c i log e -^--(c,-c,_ 1 ) 


(D-2) 


-[(8c, -5c ,)£• , - 2(5c, -2c, ,)£; +(2c i +c il )s M ]+y il +Mtan0 il 

(i = \,2,3,....,ri) 


In view of equation (D-l), equation (D-2) may be written for each index, i, as: 
For i= 1: 


[( 2 c, C () )c, C () 2(2 C, C 0 )c o £, +CjC 0 £T 2 j 

2(c,-c 0 ) 

| c i lo g e — -(c, -c 0 ) 
c o 

l-, (A/>2 [( 8 c,- 5 c,K - 2 ( 5 c, - 2 c,K +( 2 c, +, 

IZ^ c 0 J 

? oK]+yo+ A/tan6 'o 

o/ \4 [(2^2 c, )c 2 2(2 c 2 c 1 )c | c 2 +c 2 cy;. j 

2(c 2 -c,) 

C 2 l0 g« — -( C 2“ C l) 

C 1 J 

1?/ ^ ^2 [( 8c 2 — 5 c,)c; — 2 ( 5 c 2 — 2 c x )£ 2 +( 2 c 2 +c 1 )c 3 ]-l-_p 1 + A/tan#, 
IZ(C 2 

9/ [(^ C 2 C i) C 2 £ L ^(^ C 2 C l) C / ; 2 +6i 2 6 1 ^ 3 ] 

Z(^c 2 ~c x ) 

1 1 

/— \ 

<r 

1 

> w' 

1 

<r| <f 

1 1 

... .2 [(^ C 2 5 c 1 )f 1 2 ( 5 c 2 2 c, )c 2 +(2c, +c 

12(c 2 -c, ) 

iK] 


(D-3) 
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-H 


(A/) 2 


2 ( c i — c 0 ) 

(A/) 2 


[(2c, - c 0 ) c / ; <> - 2 (2c, - c 0 )C 0 /;, + c,c 0 /; 2 ] 


c,log e — -(c,-c 0 ) 


' 2(c, - c 0 ) 


[(8c, -5c 0 )c 0 -2(5c, -2c 0 )c, + (2c, +c 0 )c 2 ]+y 0 + A/tan<9 0 


(A/) 2 


2(c,-c 0 ) 3 

(A/) 2 


[(2c, - c 0 )c,c 0 - 2(2c, - c 0 )c 0 c, + c,c 0 c 2 ] log e 


4(c, c 0 ) 


[(5c, - 3c 0 )c 0 - 2(3c, - c 0 )c, + (c, + c 0 )c 2 ] + A/ tan6» 0 


A/tan^ 


After grouping terms, equation (D-4) becomes: 
(A/) 2 


y 2 




2 (c 2 -c,> 
(A/) 2 


[(2c 2 - c,)c 2 c, - 2(2 c 2 - c,)c,c 2 + c 2 c,c 3 ] 


C 2 l0 g« — ~( C 2~ C l) 


12 (c 2 -c,) 
(A/) 2 


[(8c 2 - 5c,)c, - 2(5 c 2 - 2c,)c 2 + (2c 2 + c,)c 3 ] 


-H 


2(c, ~c 0 ) 
(A/) 2 


C, log e — -(c,-C 0 ) 


i2(c,-c 0 y 
(AO 2 


-K2-1) 


[( 2c i ~ c o) c i f: o - 2(2c, -c 0 )c 0 c, +C,c 0 c 2 ] 

[(8c, - 5c 0 )/; 0 - 2(5c, - 2c 0 )c, + (2c, + c 0 )s 2 ] + y 0 
[(2c, - c 0 )c,c 0 - 2(2c, - c 0 )c 0 c, + c,c 0 c 2 ] log e 


2(c, — c 0 ) 
(A/) 2 


4<c,-c 0 y 


[(5c, -3c 0 )c 0 — 2(3c, -c 0 )c, +(c, +c 0 )c 2 ] 


0 >+(2)A/tan0 o 


For / = 3: 


y 3 = r[(2c 3 -c 2 )c 3 c 2 - 2(2c 3 -c 2 )c 2 c 3 +c 3 c 2 c 4 ] 


2(c 3 -c 2 ) 

(AO 2 


C 3l°g e — -( C 3- C 2) 


12(c 3 -c 2 f 
(A/) 2 


2(c 3 -c 2 ) 
(A/) 2 


[(8c 3 - 5c 2 )£ 2 - 2(5 c 3 - 2 c 2 )c 3 + (2c 3 + c 2 )c 4 \+y 2 +M\m0 2 


c 3 lo 8 e — -( c 3- c 2> 


12(c 3 -c 2 ) 


t[( 2c 3 - c 2> c 3^2 - 2 ( 2 c 3 - c 2 )c 2 e 3 + c 3 c 2 s 4 ] 

[(8c 3 - 5c 2 )c 2 - 2(5c 3 - 2c 2 )s 3 + (2 c 3 + c 2 )c 4 ] 


(D-4) 


(D-5) 
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( A /) 2 


2 (c 2 - Cl y 
( A /) 2 


[(2c 2 - c,)c 2 c, - 2(2c 2 - c x )c x s 2 + c 2 c x s 3 ] 


c 2 log e — ~(c 2 -c x ) 


12(c 2 -c x y 


[(8c 2 - 5 c, )/;, - 2(5 c 2 -2 c x )s 2 +(2 c 2 + c 1 )s 3 ] 


( A /) 2 


2(c,-c 0 y 

( A /) 2 

I2(c, -c o y 


[(2c, -c 0 )c,c 0 - 2(2c, -c 0 )c 0 s x +c x c 0 s 2 ] 


c,log e — -(c,-c 0 ) 


[(8c,-5c 0 )£ 0 -2(5c, -2c 0 )/;, + (2c, + c 0 )c 2 ] + _y 0 


+( 2 - 1 ) 


( A /) 2 


2(c,-c 0 ) 3 

( A /) 2 


. 4(c, -cj 
+ ( 2 )A/tan # 0 


[(2c, - c 0 )c,c 0 - 2(2c, - c 0 )c 0 £ x + c x c 0 £ 2 ] log 1 

C 0 

[(5c, - 3c 0 )c 0 - 2(3 c x - c 0 )£ x + (c, + c 0 )c 2 ] 


V2 


( A /) 2 


^(^2 ^ 1 ) 


[(2c 2 - c, )c 2 c, - 2(2c z - c x )c x £ 2 + c 2 c x £ 3 ] log e 


( A /) 2 


+4 


4(c 2 -c,y 


[(5c 2 -3c,)c, -2(3c 2 -c,)c 2 +(c 2 +c x )£ 3 ] + 


( A /) 2 


2(c, — c 0 ) 3 
( A /) 2 
4(c, -c 0 ) 2 


[(2c, - c 0 )c,c 0 - 2 (2c, - c 0 )c 0 £ x + c,c 0 c 2 ] log e ^ 

[(5c, -3c 0 )£ 0 -2(3c, -c 0 )c, + (c, + c 0 )c 2 ] + A / tan 0 o 


A/ tan 6^ 


After grouping terms, equation (D- 6 ) becomes: 

( A /) 2 


Xi = i 


2(c 3 -c 2 ) 4 
( A /) 2 


c ;i l°g fi — -(c ;s -c 2 ) 


12(c 3 -c 2 ) 
( A /) 2 


+4 


2(c 2 -c,) 
( A /) 2 


[(2c 3 - c 2 )c 3 £ 2 - 2(2 c 3 - c 2 )c 2 £ 3 + c 3 c 2 c 4 ] 

[(8c 3 - 5c 2 )c 2 - 2(5c 3 - 2c 2 )£ 3 + (2c 3 + c 2 )/; 4 ] 

c 2 l°g e — -(c 2 -c,) 


12(c 2 -c,) 


[(2c 2 - c,)c 2 c, - 2(2 c 2 - c x )c x £ 2 + c 2 c x £ 3 ] 

[(8c 2 - 5c , )£ x - 2(5 c 2 - 2 c,)c 2 + (2c 2 + c,)c 3 ] 


(D- 6 ) 
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(A /) 2 


+4 


2 (c x -c 0 y 
(A /) 2 


12 (c,-c 0 ) 

(A /) 2 


-K3-2) 


[(2c, - c 0 )c,c 0 - 2(2c, - c 0 )c 0 c, + c,c 0 c 2 ] 

[(8c, - 5c 0 )c 0 - 2(5c, - 2c 0 )c, + (2c, + c 0 )s 2 ] 
[(2c 2 - c, )c 2 c, - 2(2c 2 - c, )c, c 2 + c 2 c, ] log, 
[(5 c 2 3 c, )l\ 2(3c 2 c x )s 2 +(c 2 + c, ) c, J 


c,log e — -(c,-c 0 ) 




2(c 2 c,) 

(A /) 1 


+(3-1) 


4(c 2 c, ) 

(AO 2 


2(c, -c 0 ) 3 

(AO 2 


[(2c, -c 0 )c,£T 0 — 2(2c, -c 0 )c 0 c, + C,C 0 £ 2 ] log e — 


4(c, -cj 


-[(5c, -3c 0 )c 0 — 2(3c, -c 0 )c, +(c, +c 0 )c 2 ] 


0 [-+3A/tan0 o (D-7) 


Observing the indicial behavior, equation (D-7) can be generalized for index /, and the deflection, , 
can be expressed in dual summation form as: 


Ji-Z 

j = 1 


(A/) 2 


2(c j -c j _,> 

(A/)" 


[(2c, — Cj , XCjCj , — 2c j _,£ j ) + CjC^ ,£\ + , ] 


Cj logg “ (Cy Cy_, ) 


y-1 


12(c -c ,, ) 2 [ (8c J~ 5c J- l)g >- 1 - 2 (5c, - 2c ,, ) £j + (2c, + c,._, )c, + , ] 




Contribution from deflection terms 


+Z (*-■/) 

j = 1 


\3 [( 2c j - <y-i )0y*,-i - 2c j-x s j ) + c M jjA ] log, 
S-O yi 


4( c ( ^ t) 2 [( 5c j ~ 3 S-tH-i - 2 ( 3 S -c J -.H +(c y +c J -,)c, +1 ] 


+ (/)A/tan<9 0 


Contributions from slope terms 


(/ = 1,2, 3,...., w) (D-8) 


which is equation (27) in the text, the Improved Displacement Transfer Function for the nonuniform 
beams. 
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APPENDIX E 


DERIVATIONS OF LOG-EXPANDED SLOPE AND DEFLECTION 
EQUATIONS FOR NONUNIFORM BEAMS 

Appendix E presents the mathematical derivations of the Log-expanded slope equation (29) and the 
Log-expanded deflection equation (32). The depth factor, c(x), was considered to be a slowly-changing 
function of x. 


Log-Expanded Slope Equations 

For slowly-changing c(x) [that is,(c £ — C i l )« 1], the logarithm term, ]og e (c i i C £ l ), in the slope 
equation (22) may be expanded in series form up to the third-order term as (ref. 17): 

. c , r c , ^ i ( c , v i ( Ci v 

— = 1— - lj-~[ — -lj + ~[ — -U ----- 

<V. ) 2U m ) 3U;_, ) 


c ~ c ,-x 


C, ~ C ; 


1 g|~C,-i , (g,-Cj-i) 

2 c i \ 3c /i 

= -3c„(c, -c,_,)+2(c,-c„f 

6C i 1 

= -fe - c x-.)<3c,_.-2(c, 

6c,., 

= +(c, -c m )(2c, - 5c,.,)-....] 


(E-l) 


1-1 


In view of equation (E-l), the Improved slope equation (22) [or equation (C-17)] can be expressed as: 


tan 0, = 


M 


' 2(c,. — c f _,) 


M 

2 ( c , — c ,-,) 3 


A/ 


[(2c, - c M )(c,6' 1 , - 2c M £ i ) + cc,£ i+1 ] log„ 


“ 7-1 


~ ~ 2 C;1 [( 5 ^-3c ,V, , — 2(3c f -c, ,)c +(c, +c, Jc,,] 


-i-tan# 


[(2c, -c, ^(c-c, , - 2c l _ 1 £r-) + c,c. ,g, +1 ] M [6c* , +{c i -c, ,)(2c, - 5c, ,)] 


~ C ~ 2 il [ (5C ‘ _2 ( 3c ' + ( c < + 

+ tan0,_, 

[(2c, - c .Xcx] , - 2c i l £')+c i c i _ l £ M ][6c* l + (c, - c, , )(2c, - 5c, ,)] 

r 3c L [( 5c i ~ 3c < iK i - 2 ( 3c i - c i iH + ( c , + c , ,K„] 


12ct 1 (c,-c,_ 1 ) 2 


► + tan 0 , 
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A/ 


6c h( 2c ,- c ^) c t + c i( c i - c m)( 2c , - 5c ,-i)( 2c , - C h) 

-3cti(5c,- -3c,_,) 




4 + 


12c,_i(2c j c, j ) 2c,_,(c, c i _i)(2c I 5c,_,)(2c, c w ) 


+ tan# 


A/ 


12^ ( c , ~ c i-i) 


A1 


\2cUc-c^f 


+6c ; ,(3c i -c i _ 1 ) 

6c i c li + c i c i-i( c i ~ c «-i)( 2c i - 5 c m) 

-ScL^+c,,) 

6 + (c, -c M X4cf -12c 2 c,_, + 5c, c 2 ,) 

— 9c, ! (c, — c, | ) — 6c i _ x c i 

- ] 2c\ , 6- - 1 2cf , (c, - 6- , ) - (c, -6- ^(Scfc, , - 24c / c 2 , + 10c 3 ,) 

+12C, l c i + 6c, , (c, — c,_, ) 

+3cf_i(c I --c i _i) 

6 c ,_i c , ( c , — c ,_i ) + 6c,_,c, (c, — c,_, ) + (c, — c, | )(4c i — 1 2c, c,_, + 5c,c, , ) 
_ 9 c L( c i -C i- I ) 

+|^ — 12c,_i(c, — c,_i)— (c, — c, , )(Sc j c, , — 24c/:, , + 10c,_ 1 )+ 6c,_,(c, — c,_,)Jc, 
+[( c , - c , ,X2c - 5c, x )c,c, j + 3c 3 ,(c, - c, ,)]c, +1 


+ tan# 


+ tan# 


A/ 


12c,_ 1 (c, c,,) 


+ tan# 


A/ 


UcUc-c^) 

Al 


\} 2c L c i + 4c f ~ 12c ! c ,- l + 5c ,4i - 9c L]^-i 

+[-6<i -8c 2 c,_, +24c,c 2 , - 10c 3 ,)Jc, 

+[(2c, - + 3c 3 _,]c, +1 

[-12c,c,_i (c, - c,_, ) + 5c 2 , (c, - c,_,) + 4(c 3 - c 3 lX )]c,_, 
+[-8c,.c m (c, - c, ,)+ 16c 2 ,(c, - c, _,)]«, 

+ [ 2 (c, - c,_, )c,c, , — 3(c, - c,_, )c, 2 , ] c, +1 


+ tan 0, 


\2d 


||^ 12c, c, , + 5c, , + 4(c, + cc j + c, , )J c, , (8c, c, , 1 6c, , )c, + (2c,-c, , 3c, , )c, +1 1 


(E-2) 


+ tan# 


(/ = 1,2,3,....,k) 


After the final grouping of terms, equation (E-2) becomes: 
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taa0 i = {[ 5c L + 4 k - C /-1 ) 2 ] ^,-1 - k-i k - 2c ,-i k + C ,-1 ( 2c , - k-, k + i } + ton^-1 (E-3) 

(/ = 12,3,...., n) 

which is equation (29) in the text, the Log-expanded slope equation for the nonuniform beams. 

For the uniform beam case (c £ - = c i X — c ), equation (E-3) degenerates into: 


tan O t = ^~ (5^ + - s M , ; (i = l, 2, 3,... n) 

12 c 


(E-4) 


which is identical to equation (14) in the text. 


Log-Expanded Deflection Equations 

For the case of the Log-expanded deflection equation, the logarithm term, log (i (c i / c i ,), had to be 
expanded up to the fourth-order term as (ref. 17): 


, C , (c, P 1 ( c, V 1 ( c, V 1 ( c, V 

+ iUr‘J -jUt'j + - 

= ^= L [ 12 c?- , -64 , te-c,., ) + 4 c,. ,( c , -c,. ,)2 - 3( (c , - C ,. ,)J +....] 


(E-5) 


The term, [c f k)g < ,(c ; / c ( _, ) — (c f — C-_, )] , can then be written as: 


c,log e — -k-c,i) 

c i-l 

= C| k. — -« H )+ 4 c, -,(c, ) ! -3((c,-c,J +....]-(c, 


12c*. 

C i( C i — C i-l) 


12<, 

c ,-k ~ c n) 


12c 4 

12<. - 6 clM -c^+tc^c-c^f - 3((c i -c t J el 


12c 4 


12cf 


(Pi - C H ) - 6c ?-i k - C /-1 ) + k k - C i-1 f - 3 (k - C H ) 3 


= ( C ‘i2c* 1 ') [ UC ‘ ] ~ 6CiC ‘ l + 4C,C ‘ ' iC ‘ ~ C, l)_ 3C,(C/ ” C ‘ l)2 ] 


(E-6) 


In view of equation (E-6), equation (26) [or equation (C-21)] in the text can be rewritten as: 
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Ji 


(AO 2 


2 ^-^y 

(AQ 2 

12(c, -c t -i) 2 

(AO 2 


[(2c, - C H )(c^i - 2c I t 6; ) + c lCil s M ] 


c,log,- J —(c I -c l _ I ) 


[(8c, -5 c iX )s iX -2(5 Cj -2c ,)£• + (2c, + c, ( )^ , ] +_y, 1 + A/ tan 0_ t 


2 ( C i - c iJ 


[(c 2 +(c, -c, ,)c,)^ , - 2c, c, , +2(c, -c i _ l )c i _ x )s i +c i c i _ l e M j 
X (C n ^ l) [' 2c ; 3 _i ~ 6c, c 2 , +4c,c j l (c, -c ;1 )-3c,(c -c,,) 2 ] 


(AQ 2 

12(c,-c„) 3 

(A /) 2 

24c, 4 _j(c, -c,.,) 


12c, 

[(8c, -5c, ,)£• , - 2(5c, - 2c, , )£' + (2c, + c, , )£ n ]+y it + A/tan0, 

[( c ? +( c , -(Zcfit-i +2 (c, -c, +c,c, l£ - +1 ] 

X [l2ct, - 6c,cf_, + 4c,c, ,(c, - c,_,) - 3c, (c, - c,_,) 2 ] 
_ 2c, , [(8c, —5c, 0^', i — 2(5c,. —2c iX )£ i +(2c, +c, ,)c +1 ] 

+y,, + A/tan<9 l 


(A /) 2 


24c 4 1 (c,-c, 1 ) 2 


C i C i-I Xi 1 2c, | 6c,c, , + 4c,c, , (c, C,_, ) 3c,- (c, C, , ) 

-2c 4 t (8c, - 56,0 


/(*m) 

„3 s .. . 2 


2c ( c, , + 2(c,- - c, ,)c, , )(12c,. , - 6c,c, , + 4c,c, 1 (c, - c, x )~ 3c, (c, - c, ,) ) 
+4c i 4 ,(5c,-2c, ,) 


fO-,) 

c Pi i /; ni 42c, , — 6c, c, ! + 4c,c, 1 (c, — c, 0 — 3c,(c, — c, t ) , 
-2c 4 , (2c, + c,_,) 


/(%i) 


+ y,_i+A/tan<9, , 


0 = 1-2,3,...., n) 

(E-7) 


The calculations of the terms in equation (E-7) are shown below: 

For the f{s i X ) term: 




[ c f + ( c , — c ,-i ) c i ) ( l2c h - 6c ^i + 4c,c, \i c i — c , i ) - 3c /( c / — c n f ) 

-2^(8^ -5c m ) 
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i2cfei ■ -t&L + *Zc„(c t 

+12c i c/_ 1 (c J - — Cj.j) — 6 cfc*j(c f -c,0 + 4c 2 c, ] (c, -c H ) 2 -3c 2 (c, -c^,) 3 
-lecc^ + Kkf, 

~~ ^ C i C i-l( C i _ C /-l ) + ^ C i C !-l( C i — C i- 1) _ 3c f ( C i — C i-\ ) 

+\2cfl x {c l - c,, ) - 6c 2 c 2 , (c, - c, x ) + 4c 2 c,, (c, -c, ,) 2 -3c 2 (c,- -c M ) 3 

+ £ >C 1 C ] \ ( c , - C ,1 ) - 1 ° C , 4 1 ( C z - C /-l) 

-6cfc 2 ! + 4c 3 c,_, - 3c 3 (c, -c H ) 

+1 2c, c, j — 6c, Cj_j + 4c c ! (c, — c f _, ) — 3c, (c, — c f _, ) x 




'/I 


= ( c ,- c / 1) 


= ( c ,- “<?-,) 


-t^.-lOc 4 , 

2 / \ , yl 3 o 3/ 


1 2cc il (c i C j ) + 4c c , 3c (c, Cj_j) 
44cfc M ( Cl - -c H ) - 3c, 2 (c,-c, j) 2 
“^cc 3 , + 1 Occ 3 , — 1 0c, 4 , 


^ i 


12cc, ,(c, c, , ) + 4c, c, , (c, c,_j) 3c,. (c, c, ,) 

= ( c i~ c t-i) +^ c ] c i-i( c t ~ Cj_j ) — 3c, 2 (c, — c 4 _j ) 2 
+10c 3 1 (c-c,_ ] ) 

= (c, — c, j ) |^— 12^^,^ + 4c,c, j (c, + c, j ) — 3c, + 4c, c, x — 3c, (c, — c,_ x ) + 1 0 c, x J c, x 
= (c f — c Vl ) [7c^j + 8c,c^j(c, — c,.,) — 3c,. (c,. — c M ) — 3(c, — j 

After grouping terms, one obtains the final form of f (/; , ) as: 

/<*,-, ) = (<i - C„, ) 2 [ 7<£, + (c, - c,_, XSc.c,., - 3c, 2 ) - 3(c- - c ’_, )] s M 


For the /(^r ) term: 


/(«i)= 


-{2c,-c,_i +2(c, -c,-_ 1 )c,)(l2c 3 _ 1 -6c, +4c,c,_ 1 (c, -c,_ 1 )-3c,(c,- -c,_,) 2 ) 
— 2c, 4 _j (2(5c, — 2c, _ x )) 

-24c, c, 4 , — 24c 4 (c, — Cj_j ) + 1 2c, 2 c 3 + 12c, c, 3 (c, -c M ) ' 

-8c 2 c 2 i(c, - c, j ) - 8c,r, 2 j (c, - c (1 ) 2 
+6c,c, j(c, — c, j) + 6c i c f _ 1 (c, — c f _j) + 20c,c, j — 8c, j 


-^c 4 , +12^-8*^ 


( c , ~Ci- 1) 


-24c 4 + 12c, c, 3 -8 c)c]_ x -c, t ) 

+6c 2 c i l (c - c, ,) + 6c, c,. ,(c, -c,. t ) 2 


>-c. 


(E-8) 


(E-9) 
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-c,-_i)+8t£i(cf -c 2 ,) 

-24c, 1 + 12c,c 3 -8c, 2 c 2 - 8c,c, 2 , (c, -c,,) 
+6c ; 2 c i _ 1 (c, - c f _, ) + 6c,c,_ 1 (c, - c H ) 2 
4c i c^ 1 + 8c 3 1 (c, + c, 1 )-24< 1 + 12^, -8c 2 c 2 , 
-8c,4i ( c , - c «-i ) + 6c f c ,-i ( c , ~ <Vi ) + 6^,. (Pi ~ c,-x f 


( c , 

= ( c , - c /-i) 
= ( c i- c ,-i) 
= ( c , - c /-i) 


1 6c, c, 16c,! + 8c,c. , 8c, c, 

— 8c,c, , (c, — c, , ) + 6c, c, j (c, — c, ) + 6c,c, , (c, — c,_ t ) 


1 6c, i (c,. - c, j ) - 8c,.c ,- 1 (c, - c, , ) 

-8 cfl, (c, - c,_! ) + 6c 2 c,_, (c, - c,, t ) + 6c, c, , (c, - c,_, ) 2 
= (c, — c, j) p6c, , — 16c,c, j + 6c, c, , + 6c,c, , (c, — c, j)Jc, 

= (c, -C, ,) 2 (16c 3 ] - 22c, c 2 +12c 2 c, ,)c, 

= 2c, ,(c, - c, ]) 2 (6c, 2 - 1 lc,c, , + 8c 2 ,)c, 


After grouping terms, one obtains the following final form of f(s t ) as: 

/(«,) = 2c,., (c, - c,_, ) 2 [3c? + (c, - c,., )(3c, - 8c,_, )>, 


For the /(c, fl ) term: 


/(%,)- 


c,c, i 12c, , - 6c, c, , + 4c, c, ,(c, -c, ,)- 3c,(c,-c, ,) ) 
_-2c I 4 _ 1 (2c, +c,._j) 

12c, c, 4 ., - 6c 2 c 3 + 4c 2 c 2 1 (c,. - c, , ) - 3c 2 c, t (c f - c, , ) 2 


"7+1 




2c 


i-l 

2 3 r> 5 




7+1 


> £ m 


= 8c,< 1 -6c 2 c;-2c^ 1 

+(c, - c i _,)[4c 2 c 2 j - 3c 2 c,_ 1 (c, - c,,)] 

-6c i cti(c i -c i _ 1 ) + 2c 4 I (c i -c j _ 1 ) 

+(c, - c i _,)[4c 2 c 2 j - 3c 2 c,_ 1 (c, - c,,)] 

= (c, — c, i)j — 6c,c’, , + 2c, j + 4c, c,_, — 3c, c, , (c, — c, , ) jc J+1 
= (c, — c, i ) [^ 2c, , (c, — c, , ) + 4c, c, , (c, — c, ,) — 3c, c, ,(c, — c, , )Jc, 
= ( c / -c,_i) 2 (-2c 3 _, + 4c, c 2 , - 3c 2 c, ,)c, +1 


After grouping terms, equation (E-12) becomes: 

/(% 1 ) = - c i-i( c / - c i-i) 2 [ c ? + 2(c, - c, ,) 2 ]c,_, 


(E-10) 


(E-ll) 


(E-12) 


(E-13) 
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In view of equations (E-9), (E-l 1), and (El 3), the displacement equation (E-7) becomes: 

(A/) 2 


* 24 cUc-c^f 


( A /) 2 

2AcU Cl - Cl J 


1 ) + 1 ) + /(*i_i )] + y t . 1 + V tan^ j 

( c , - c , i) 2 [ 7c L + ( c , - c i 1 X 8 c fi I - 3c 2 ) - 3(cf - c ] , )]*• , 

< +2c i l (c i -c i l ) 2 [(3c 2 +(c i -c i l )(3c i -8c i ,)]^. 

— C i-l( C i — C i-l) ^C,- + 2(6^ — C ; , ) Jc j+1 

+_y i( + A/ tan , 

(/ = 1,2,3,....,/?) 


(E-14) 


After eliminating the factor (c, -C, ,) in the denominator, there results: 

_ (A/) 2 J [ 7c L + ( c i ~ c.-i X^A‘, i -3c 2 )- 3(c* - )]*•,. , 

24c/_ t +2Cii [(3c 2 + (c, - X3c, - 8c M )]*, - c M [c 2 + 2(c f -c M ) 2 ]* l4 
+y ( _, + A/tant^, 

(/ = 1,2,3,....,/?) 

which is equation (32) in the text, the Log-expanded deflection equation for the nonuniform beams. 


(E-l 5) 
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APPENDIX F 


DERIVATIONS OF LOG-EXPANDED DISPLACEMENT TRANSFER 
FUNCTION FOR NONUNIFORM BEAMS 

Appendix F presents the derivation of the final Log-expanded Displacement Transfer Function 
written in dual summation form. The Log-expanded slope equation (29) [or equation (E-3)] and the 
Log-expanded deflection equation (32) [or equation (E-15)] are duplicated below as equations (F-l) and 
(F-2), respectively: 


tan 4 = {[ 5c t i + 4<c, - - 8c f _i(c t - 2c, ,)/.• + c- ,(2c, - 3 c, ,)c i+1 } + tan 6>_, (F-l) 

(* = U3,...,n) 


y, 


( A /) 2 
24c, 4 , 


[7c^ + (c, - - 3c 2 ) - 3(c 3 

_ + 2c ,-i [(3c 2 + (c - c , )(3c - Sc; , )] c - c,_, [c 2 +2(c i -c i _ l f]e^ 

+y i _ l + A/ tan 0 t , 

0 = 1,2,3,....,«) 


(F-2) 


In view of equation (F-l), equation (F-2) can be written out explicitly for different values of the 
index, /, as follows. 


For / = 1 : 


(A /) 2 CoX&W) O>o 

24c 4 +2c 0 [(3c 1 2 + (c, - c 0 )(3c t - 8c 0 )] - c 0 [c 2 + 2(c t -c 0 ) 2 ]c 2 

+y 0 + A/tan # 0 


(F-3) 


For z = 2: 


y 2 


( A /) 2 

24c 4 


[7c, + (c 2 - c,)(8c 2 c, - 3c 2 ) - 3(c 2 - c 3 )]c, 

+2Cj [(3c 2 + (c 2 - Cj)(3c 2 - 8c,)]c 2 - c, [c 2 + 2(c 2 - c,) 2 ]c 3 
+ y, + A/tan0, 


Substitutions of equations (F-l) (for Z = 1) and (F-3) into equation (F-4) yield: 
(A /) 2 


y 2 = 


24c, 4 


[7c, 3 + (c 2 - c, X8 c 2 c, - 3c 2 ) - 3 (c 2 - c 3 )]c, 

+2c,[(3c 2 +(c 2 -c,)(3c 2 -8c,)]^ 2 -c,[c 2 + 2(c 2 -c,) 2 ]c 3 


(F-4) 
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>+j> 0 +A/tan0 o 


+ 


(A If 

24 ci 


[?c o +(c, -c 0 )(8c,c 0 -3cj 2 )-3(cj 3 -c 3 )]^ 0 
+2c 0 [(3cj + (c, — c 0 )(3c, — 8c 0 )J£, — c 0 [c, + 2(c,— c 0 ) J c 


(A/) 2 

12c 3 


{[5c 0 2 +4(c, -c 0 ) 2 ]c 0 -8c 0 (c, -2c 0 )c, + c 0 (2c, - 3c 0 )c 2 } + At tan 0 0 


(F-5) 


Ai t;in (\ 


y 2 = 


After grouping terms, equation (F-5) becomes: 

(A/) 


24c, 4 

(A/) 2 

24c 4 


\lcl + (c 2 -c x )(Sc 2 c x -3c 2 2 )-3(c 2 3 

|^(3c 2 + (c 2 — )(3 c 2 — ) J^ 2 — c x | ~ c 2 + 2(c 2 — ) J s 3 

Oo + (^*i — X&W, — -^i ) — — c 0 )J s Q 

r* 

+2c 0 + ( c x — c 0 )(3 c x — 8c 0 ) J s x — c Q | ~c x + 2 (c x — c 0 ) J^ 2 

+ (2 - {[5c 2 + 4(c, - c 0 ) 2 ]c 0 - 8c 0 (c, - 2c 0 )c, + c 0 (2c, - 3c 0 )c 2 } + (2)A/tan0 o 

Yz,C n 


+y 0 


(F-6) 


For i = 3: 


y 3 = 


(A If 
24 c 4 


[7c 3 +(c 3 -c 2 )(8c 3 c 2 -3c 2 )-3(c 3 -c 3 )]c 2 
+2c 2 [(3c 2 +(c 3 -c 2 )(3c 3 -8c 2 )]c 3 -c 2 [c 2 +2(c 3 -c 2 ) 2 ]c 4 
+y 2 + Altaa0 2 


(F-7) 


Substitutions of equations (F-l) (for i = 2) and (F-6) into equation (F-7) yield: 

(A/) 2 


y 3 


24c 4 


(A/) 2 

24c 4 


[7c 3 +(c 3 — c 2 )(8 c 3 c 2 — 3c 2 ) — 3(c 3 — c 3 )] s 2 

+2c 2 [(3c. + (c 3 — c 2 )(3c 3 8c, )J c. — c 2 [c 3 +2(c 3 — c 2 ) Jc 

[7c 3 + (c 2 - c, )(8 c 2 c, - 3c 2 ) - 3(c 3 - c, )] c, 

+2c, [(3c 2 + (c 2 - c,)(3c 2 - 8c,)]c 2 - c, [c 2 + 2(c 2 - c,) 2 ]c 3 


+ -Fo 


(A If pci +(c, -c 0 )(8c,c 0 -3c 2 )-3(c 3 -c 3 )]c 0 
24c 4 1+2 c 0 [(3c, 2 + (c, - c 0 )(3c, - 8c 0 )]c, - c 0 [c 2 + 2(c, - c 0 ) 2 ] c 2 

+(2-l)^V{[ 5c o + 4 ( c i -c 0 f^\s 0 ~ 8c 0 (c, — 2c 0 )c, + c 0 (2c, - 3c 0 )c 2 } + (2)A/ tan 0 a 


y 2 
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+ ^l2c 3 {[ 5C,2+4 ^ 2 _C l) 2 ]^ _8C l( C 2 - 2C lK + C i( 2C 2 - 3C lH} 

+ ^r{[ 5c o +4( C 1 -<?o) 2 ]^o -8c 0 ( c i - 2c oK + c o( 2c i -3 c 0 K} + A/tan6> 0 


(F-8) 


AI t;in 0 ? 


After grouping terms, equation (F-8) becomes: 

(A/) 2 


y 3 


24c 2 


(A If 

24c! 


(A/) 2 

24c! 


\j c 2 + ( c 3 c 2 )(8c 3 c 2 3 c 3 ) 3(c 3 c 2 )J^ 2 

+2c 2 [(3c 2 + (c 3 - c 2 )(3c 3 - 8 c 2 )]c 3 - c 2 [c 2 + 2 (c 3 - c 2 ) 2 ]c 4 
[ 2c i 3" (c 2 — Cj )(8c 2 c, — 3c 2 ) — 3(c 2 — c, ) J c, 

+ 2 c \ [(3c 2 + (c 2 -c,)(3c 2 -8c,)]c 2 -c, [c 2 + 2(c 2 -c,) 2 ]s 3 
[7c 0 + (c, — c 0 )(8CjC 0 — 3Cj ) — 3(c, c Q ) J £ 0 
+2c o [( 3c i 3"( c i — c 0 )(3c, — 8c 0 )Jc, — c 0 [c, 3" 2(c, — C„) J *'i 

+( 3 - 2 )^r {[ 5c i 2 + 4 ( c 2 - C 1 ) 2 ] /; ! - 8c, (c 2 - 2c, )/; 2 + c, (2c 2 - 3c, ) s 3 } 
3-( 3 - 1 )^t{[ 5c o + 4 ( c i -„) 2 > 0 -8c 0 (c, -2c 0 )c, +c 0 (2c, - 3c 0 )c 2 } + (3)A/ tan6> 0 


3- To 


(F-9) 


Observing the indicial behavior, equation (F-9) can be generalized for index /, and the deflection, _y 2 , 
can be expressed in dual summation form as: 


(A Iff 1 

24 M C U 

[7c’_, _,)(8c jCj _, -3c])-3(c] 

+ 2c /-i [( 3c y + ( c j _ e ,-iX 3c y “ Scj-i )]«; - «j-i [<v 

-1 

+ 2 (Cj ~ Cj j) 

>+y 0 


V 

Contributions from deflection terms 



(A/) 2 y, (i-j) f[ 5c /-i 3“ 4 (c, Cj i) ]c,_, 8c y , (Cj 2c J _ l )s J 
12 M c j - 1 1 +c j _ l (2c j -3c J _ l )e j+l 

> + (j)A/ tan 0 O 



Contributions from slope terms 

(/ = 1,2,3,...., n) (F-10) 


which is equation (33) in the text, the Log-expanded Displacement Transfer Function for the nonuniform 
beams. 
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APPENDIX G 


DERIVATIONS OF DEPTH-EXPANDED SLOPE AND DEFLECTION 
EQUATIONS FOR NONUNIFORM BEAMS 

Appendix G presents the details of the derivations of the Depth-expanded slope equation (41) and the 
Depth-expanded deflection equation (46) for nonuniform beams. The reciprocal of the depth factor, 
l/c(x), was first expanded in series form before carrying out the integrations of the curvature 
equation (3) in the text. 


Depth-Expanded Slope Equation 

The slope, tan#(x)D of the nonuniform beam in the region, x tl < x < x i ? between the two adjacent 
strain-sensing stations, { X i l ,X i }, can be obtained by integrating the nonuniform-beam curvature 

equation (3) with the constant of integration determined by enforcing the continuity of slope at the 
inboard adjacent strain-sensing station, x 2l , as 

ia " ^ = c 57* + —fa = r a + tan (G-i) 

v v J Slope at Xj j » * > 

Slope increment gq (3) 

The distributions of beam depth factor, c{x) , and bending strain, £(x), in the domain, x i X < x < x /? 
between the two adjacent strain-sensing stations, { X i l ,X i }, described by equation (10) are rewritten 
respectively below: 


c(x) = c, , + (c, - c H ) * ^ ; (x,, < x < x, ) 

A/ 


(G-2) 


1 1 2 

£(x) = £ tl - —(3^, , - 4s i + s M )(x - x M ) + ^ (,£, , - 2s t + e M )(x - X, ,) 


(x ;1 <x<x,) (G-3) 


Substituting equations (G-2) and (G-3) into equation (G-l) yields: 


tan 0(x) = f 

Jx 


5 i - (3 £)-i - 4c, + c M ) + (e H -2c, + s l+1 ) 


c i \ +(g ~ c i .) 


x — x 

2 A/ 


dx + tan 61, (G-4) 

(x, , <x<x,) 


Because the slope term, (c i —c i l )/Al, of the depth factor, c(x), in equation (G-2) is small, its 
reciprocal, 1 / c(x), can be expanded in binomial series form as (ref. 17): 
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1 

c(jc) 


1 


1 


c ,-, +(^-0" A/ ,_i 


i+ 


( c, t x - : 


ur'j 


A/ 


1 


, ( c .Ix-r. ( c V (x-x f 

1 — I — - — 1 1 — + 1 — — 1| l -=r- 

\ c t _ x ) A/ \ c tl ) (A/) 


(G-5) 


<X-, <x<x,) 


Let 


A =_^t-l~ 4S i +S M 


B 


2M 


_s l _ l -2s l +s H 


2(A/) 2 


(G-6) 

(G-7) 

(G-8) 

(G-9) 


J A/ 

4 = (x-x, t ) 

then, equation (G-4) takes on the following simplified form: 

mn «(#) = — fe.+Jf+BfXt-Df+jyf- ....)cl4+*m0,_, 

Gi Jo 

= — ! -e^DZ+e^I?? +A£-AD£ 2 +AD 1 ? +B? -BDf +....]</£+ tan 0, , 

C l 1 ° 

= — -(A- £ t _ x D)^ + {-AD + B + s iX D 2 )g 2 - (BD - AD 2 )g 3 +....]</£ + tan 6> t 

C i - 1 ° 

(0<£<A/) (G-10) 

After carrying out the integration of equation (G-10), one obtains (ref. 17): 

1 


tan0(4 r ) = 


z-l 


s,_ l 4-(A-e,_fi)£-+(-AD+B+s hl D‘)^-(BD-AD‘)^+.... 


+ lan ^ 


z- 1 


(0<£<A I) (G- 11) 


In view of the definitions from equations (G-6), (G-7), (G-8), and (G-9), equation (G-l 1) takes on the 
following form: 
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tan<9(£)= ' \s t ^~ 


(3c -4c+c ) ( c 1 

- 5=1 ' —+£; ,1 — “I I 

2A1 ' 'U,, )M 


I 2 

2 


2M 




A / 


2(A If 


i Ml. £ \ I 

2(Uf U,_, ) 


A/ 


c, ,Y 

i 

c,_, ) (A/) 2 

\ 2 

l) - ‘ 

<r > 

J m 

4 


3 

-i-tan# 


(0<£<A/) (G-12) 

At the strain-sensing station x r one can write £, (x, x, ,) A/, and equation (G-12) gives the 

slope, tan 0 i [= tan#(x,)] , at the strain-sensing station, x t as: 


tan# ; = 


M 


L„-i 

(3c -4c Ac) f c 

[ 11 2 

2 A/ ‘“V M J_ 


1 

+ - 
3 


(3e H -4e,+e m )( c, j (i\ , -2s t +£;„) ( c, V 

2 U,, J 2 He,, J 


(Sj_i 2 c, + c,, ,) ( c, _ j'j _j_ (3fi-_i 4^-+c , ri ) ( c, _ j'j 

2 U M J 2A l U M J 


>+tan# 


A/ 

12c~ 


!l2c,_,-6 


(3s i l -4s i +s M ) ( c, ,"l 
2 


+ 4 


-3 


(3^,-4^ + ^jf c, "j (s_ l -2s i + s i J ( C' V 

2 W, J 2 '"U, J 


(s, i ~2s l +£ l+l ) ( c t y | (3s i X -4s i + s M ) f _c L _ 1 V 

2 U,_ t J 2 U, , J 


► + tan# 


A/ 

12c. 


{l2^-i - 3 (3c, , - 4c, + £ m ) + 2(^Vi - 2c, + c, +1 ) 


r Ci ^ 

+ ur‘J 

A, ,V 

+ ur‘j 


-6c, + 2(3c,_, - 4c, + e M ) - - (c,_, - 2c, + c, +1 ) 


4 ^-i--(3^-,-4c, + c, +1 ) 


+ tan 0, 


66 


(G-13) 


A/ 

I2c~ 


(5 £ i-i + % £ i ~ £ M ) + ~ 1 j { ~ l £ i-i ~ $ £ i + l £ , 


( c t Vf 1 , 3 ) 

Ur 1 ) l-2^ +6 *-2**) 


+tan# 


(i = U2, 3,. 


After rearrangements, equation (G-13) becomes: 

A/ 


tan 6, = 


' 24 


( c } 

2(5^ + Se t - s i+x ) - (-*- - lj (3^ +10^ - s M ) 


( Cl V 




(G-14) 


+tan# 


-1| {s l _ l -\l£ i +'i£ M ) 

(* = 12, 3,. 

which is equation (41) in the text, the Depth-expanded slope equation for the nonuniform beams. 


Depth-Expanded Deflection Equations 

The deflection, y(x ), of the uniform beam in the small domain, x i X <x<x L , between the two 
adjacent strain-sensing stations, { can be obtained by integrating the slope equation (G-12) with 

constant of integration determined by enforcing the continuity of deflection at the inboard adjacent strain- 
sensing station, x t l , as 


X^) = f X tan^r)rA+ = ^\znO(^)dZ + y t , ; (x ;1 <x<x t ) 


Eq. (G-12) 


Deflection Eq. (G-12) 


(G-15) 


Substituting equation (G-12) into equation (G-15), and carrying out integrations, one obtains (ref. 17): 


y(4) = — i g M — - 


(3 g ,-i -4g i + g i + i) , „ l 

2 A / ‘"V H )M 


£ 

6 


+ 


( 3 5-i ~ + £ m ) ( c Y 1 | ( £ i -2 g f + £ lt ) | g jV 1 


2AI 


xtr) 


A / 


2(A/) 


Ur 1 ) <*? 


£_ 

12 


(3 g / 1 — 4 £ , + e M ) ( c- 

l^j 1 _,_(^-i + ^ /+1 ) ( c i 

,) 1 


2 A/ U_ t 

'J (AO 2 2(A/) 2 U„ 

A/ 

20 


+3 ; ,i+^tan6( 1 


(0<^< A/) (G-16) 
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At the strain-sensing station JC i? one can write g = x i —x i = A/, and equation (G-16) yields the 
deflection, y z [= X*,)]? at the strain-sensing station, x /5 as follows: 


(A/) 2 1 1 

^ 772^-6 


(3*,— i— ' 4$+$ +1 ) , „ f $ A 

2 '-‘U,_rJ 


12 

1 

20 


(3^1 + £ i+l ) f c i {£ t j 2/; + £: /+1 ) f c z il 

2 urJ 2 + *-urJ 


(3^_ 1 -4g t +fiJ [ r c, A | ($_,- 2c, + c, +1 ) t d 

2 U,_ t ) 2 U M J 


+y ;1 + A/tan6>, 


(AO 2 


60c, 


30£ M -5(3 £M -4£, + £ M )-10^-lj« H 


5 

2 ' 


f c, t 5 


/cl .V 


+ -(•^-1 -4$ + +-($-i -2$ + $J + 5[^-- lj $_, 


tc. .V 3 


(c, A 


- ^(3$_i —4^ + ^)^—= — lj --($ t -2$ + $ +1 )| —-I 




+y [X + A/tan6( t 


(A /) 2 


60c 


30c ;1 - 5(3c,_, - 4c, + $ +1 )+-(c,, - 2c, + c, +1 ) 


+ \J*~ ~ lj - 10 $t + 1 (3 $_, - 4c, + c, +1 >- |($_, - 2c, + c w ) 
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After rearranging equation (G-17), there results: 


y-,= 


(Alf 


24c, 


2 ( c } 

0 £ i~\ +es i -s M )--y^—\^4£ i _ l + l£-s M ) 




< 0 - 


+y ( _, + A/tan O j , 

0 = 1, 2, 3,. 


which is equation (46) in the text, the Depth-expanded deflection equation for the nonuniform beams. 


APPENDIX H 


DERIVATIONS OF DEPTH-EXPANDED DISPLACEMENT TRANSFER 
FUNCTION FOR NONUNIFORM BEAMS 

Appendix H presents the derivation of the final Depth-expanded Displacement Transfer Function 
written in dual summation form. The Depth-expanded slope equation (41) [or equation (G-14)] and 
Depth-expanded deflection equation (46) [or equation (G-18)] are duplicated below as equations (H-l) 
and (H-2) respectively: 


tan 0. = 


Al 


‘ 24c ; 


2(5e„ +fe, l)(fe„ +10£, <%. 


+ tan# 


0 = (H-l) 


y,= 


(Al ) 2 


24 c. 


2 fa ^ 


if a V 


(7 £ i t +6 £ i 1J <4^/-! + 7 ^-^„) + 5 ^ c ' -lj («i_i +I2e,-3£ M ) 


+y t ! + Allan 0 iX 


(/ = L2,3, ---,«) (H-2) 


In view of equation (H-l), equation (H-2) may be written out explicitly for different values of index, 
i, as follows: 


For / = 1 : 
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(H-3) 


For / = 2: 
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+< 
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2 ( c } ifc V 

(ls 0 + 6c, - c 2 ) - - lj (4c 0 + 7c, - c 2 ) + -(^ - lj (c 0 +12c, - 3c 2 ) 

+_y 0 + A/tan6> 0 


(A/) 2 


24c„ 


2(5c 0 + 8c, - c 2 ) - ^ - lj (3c 0 + 10c, - c 2 ) - ^ - lj (c 0 - 12c, + 3c 2 ) 


+ A/tan#„ 


AI tan#, 


(H-4) 


After grouping terms equation (H-4) becomes: 
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For z = 3: 
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+4 
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After grouping terms, equation (H-6) becomes: 
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Observing the indicial behavior, equation (H-7) can be generalized for index i, and the deflection, y t , 
can be expressed in dual summation form as: 
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Contributions from deflection terms 
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Contributions from slope te rmss 


0 = 1,2,3,...., w) 


which is equation (48) in the text, the Depth-expanded Displacement Transfer Function for the 
nonuniform beams. 
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APPENDIX I 


FLOW CHART FOR STRUCTURE DEFORMED SHAPE 

VISUALIZATIONS 

(cf., U.S. Patent No. 7,520,176, issued April, 2009) 
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(Two-line strain-sensing system, refs. 1-3). 


Displacement Transfer Functions (Depth-Expanded): 


For the front strain-sensing line: 
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For the rear strain-sensing line: 
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Contributions from slope terms 


0 = U2,3,....,k) 


Cross-Sectional Twist Angle Equation: 

^sin'j^^j . (/ = 0,1,2,3,....,») (1-3) 

(4= strain-sensing lines separation distance, fig. 1-1) 

Tapered Wing Box 

Figure 1-1 shows the application of the two-line strain-sensing system to a tapered wing box structure 
(the neutral surface is located at half depth) to sense the surface bending strains induced by combined 
bending and torsion loadings. 
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c i = c O-( c O- c r,)ji 



Figure 1-1. Tapered cantilever wing box instrumented with two-line strain-sensing system for deformed 
shape predictions of structure under combined bending and torsion (ref. 2). 

The two strain-sensing lines (front and rear) are located along the lower edges (front and rear) of the 
wing box with chord- wise separation distance, d v which decreases toward the wing box tip according to 

the wing box horizontal taper rate. The two strain-sensing lines may also be aligned parallel to each other 
(that is, d= constant) along the span of the wing box. 

As shown in the flow chart, the surface strains, = 0, \ 2, 3,...., ri) are to be input to the left box of 

the chart, and the Displacement Transfer Functions {equation (1-1), equation (1-2)} (the second box from 
the left) are to convert the surface strains into out-of-plane deflections, {y r y '} (the third box from the 

left), from which the cross-sectional twist angles, can be calculated from the cross-sectional twist 

angle equation, equation (1-3). The resulting data could then be used to map the overall structural 
deformed shapes for visual display using the deformed-shape-visualization computer program (the fourth 
box from the left), and for control feedback. 

If the neutral surface is not at the half depth, a second two-line strain-sensing system must be installed 
on the upper surface, right on top of the lower sensing lines. Using pairs of lower and upper strains, the 
depth factor, c.(i = 0,l,2,3,.... ? /2), needed for the deflection calculations, can be determined (see the 

Ikhana wing shape predictions, ref. 3). 

Free-Free Beam Structure 

A good example of the free-free (unsupported) beam structure is an aircraft fuselage during flight 
(fig. 1-2). 
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Strain sensor 


Horizontal bending 
strain-sensing line 



Figure 1-2. Two-line strain-sensing system installed on an aircraft fuselage for horizontal and vertical 
bending strain sensing (ref 2). 

The fuselage may be considered as an union of two nonuniform cantilever tubular beams (that is, a 
fore- and an aft-fuselage section) that are joined together at the center-of-gravity (CG) cross section 
(fig. 1-2). For monitoring the deflections of the airborne fuselage under vertical bending and horizontal 
bending, two strain-sensing lines can be installed axially on the bottom (belly) surface and on the side of 
the fuselage (fig. 1-2). During the flight, the CG cross section has the least movement compared with the 
other fuselage cross sections. For practical purposes, it is best to choose the reference strain-sensing 
station, x 0 , at the CG cross section, and divide each strain-sensing line into two segments (fig. 1-2). One 

segment is for the fore-fuselage section and the other for the aft-fuselage section. Thus, by imposing 
y 0 =tan # 0 =0 at the reference sensing station x 0 (considered as a moving built-in end), deflection 

equation (1-1) can be used for the bottom strain-sensing line to calculate the vertical deflections of the 
fore- and aft-fuselage sections. The vertical deflection, y-, calculated from equation (1-1) will then be the 

relative deflection with respect to the tangent line passing though the reference sensing station, x 0 , which 

moves with the fuselage during flight (fig. 1-2). The same arguments hold for the strain-sensing line for 
the horizontal bending. 
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